﻿TEOREMA DIVERGENȚEI DACĂ • n dS = Ш, V-FdV TEOREMA LUI STOKES F • t ds s IDENTĂȚI ÎN IMPLICAȚI OPERATORUL V* V(/g) =fVg + gVf V(F • G) = (G • V)F + (F • V)G + FX (VXG) + GX (VXF) V-(/F)=/V- + FV/' V • (FXG) = G • (VXF) - F • (VXG) V VXF = VX (/F) =/VXF + (V/) XF VX (FXG) = (G • V)F - (F • V)G + F(V • G) - G(V • F) VX (VXF) = V(V • F) - V F VXV/ = /și g sunt funcții scalare de poziție, iar F și G sunt funcții vectoriale de poziție div, grad, curl și toată a patra ediție DIVISIUNEA DE AFACERI/ȘTIINȚĂ/TEHNOLOGIE BIBLIOTECA PUBLICĂ CHICAGO SOUTH STATE STREET CHICAGO, IL div, absolvent, un text informativ hm schey răsuci, si toate astea pe calculul vectorial a patra editie W • W • NORTON & COMPANY New York • Londra WW Norton & Company este independentă de la înființarea sa în , când William Warder Norton și Mary D Herter Norton au publicat pentru prima dată prelegeri susținute la Institutul Poporului, divizia de educație a adulților a Cooper Union din New York City Soții Norton și-au extins în curând programul dincolo de Institut, publicând cărți ale unor academicieni celebri din America și din străinătate Până la mijlocul secolului, cei doi piloni majori ai programului de publicare al lui Norton - cărțile de comerț și textele pentru colegii - erau ferm stabiliți În anii , familia Norton a transferat controlul companiei către angajații săi, iar astăzi, cu un personal de patru sute și un număr comparabil de titluri profesionale, universitare și profesionale publicate în fiecare an, WW Norton & Company este cea mai mare și cea mai veche editura detinuta in totalitate de angajatii sai Copyright (c) , , de către WW Norton & Compania, Inc Toate drepturile rezervate Tipărit în Statele Unite ale Americii Compoziție de GGS Book Services, Atlantic Highlands Datele de catalogare în publicație ale Bibliotecii Congresului Schey, HM (Harry Moritz), - Div, grad, curl, and all that: an information text on vector calculus / HM Schey -ed a IV-a p cm Include index Analiza vectoriala I Titlu QA S ' -dc ISBN - - - (pbk ) WW Norton & Company Inc , Fifth Avenue, New York, NY www wwnorton com WW Norton & Company Ltd , Castle House, / Wells Street, Londra W T QT RD D ччаза SUSWESS/SCITNCE/TECHNC • MINCIUNEA PUBLICA CHICAGO SOUTH STATE STREET снюлео, jl боеоб Zusammengestohlen aus Verschiedenem diesem und jenem Ludwig van Beethoven Cuprins Prefață la ediția a patra ix Capitolul I Introducere, Funcții vectoriale și Electrostatică Introducere Funcții vectoriale Electrostatic Probleme Capitolul II Integrale de suprafață și Divergenta Legea lui Gauss Vectorul normal al unității Definiția integralelor de suprafață Evaluarea integralelor de suprafață Fluxul Folosirea legii lui Gauss pentru a găsi câmpul Divergenta Divergența între cilindric și sferic Coordonatele Notația Del Teorema divergenței Două aplicații simple ale divergenței Teorema Probleme Contente Capitolul III Integrale de linie și bucla Integrale de lucru și linie Integrale de linie care implică funcții vectoriale Calea Independentei The Curl Curl în coordonate cilindrice și sferice Semnificația buclei Forma diferențiată a Legii circulației Teorema lui Stokes ' O aplicare a teoremei lui Stokes Teorema lui Stokes și regiunile pur și simplu conectate Independența căii și bucla Probleme Capitolul ГѴ Gradientul Integrale de linie și gradientul Găsirea câmpului electrostatic Folosind ecuația lui Laplace Direcțional Derivatives and the Gradient Semnificația geometrică a gradientului Gradientul în cilindric și sferic Coordonatele Probleme Soluții la probleme Index viii Prefață la ediția a patra Putem avea vreodată prea multe lucruri bune? Miguel de Cervantes Această nouă ediție diferă de a treia în două aspecte majore În primul rând, au fost adăugate o serie de exemple noi lucrate Acest lucru a fost făcut ca răspuns la comentariile multor studenți că astfel de exemple ar fi un ajutor în înțelegerea materialului și ar fi utile în pregătirea lor pentru a face seturile de probleme Sarcina mea a fost să adaug suficiente exemple pentru a fi de ajutor, fără a prelungi în același timp textul în mod semnificativ (Doi recenzenți m-au îndemnat să nu prelungesc deloc cartea, oferind din neatenție un răspuns la întrebarea Sr de Cervantes ) A doua diferență majoră între această ediție și predecesorul ei implică schimbarea rolurilor celor două unghiuri sferice Ѳ și ф Cu ani în urmă, când a fost scrisă cartea, era obișnuit să se folosească Ѳ ca unghi polar și ф ca azimutal În zilele noastre, convenția mai comună inversează acest lucru, făcând Ѳ unghiul azimutal și ф polar, convenția pe care o adoptăm în această ediție Doresc să le mulțumesc numeroșilor cititori care, de-a lungul anilor, mi-au scris cu sugestii de îmbunătățire a textului Aceste sugestii au fost adesea adoptate și au fost un motiv important pentru durata foarte lungă de viață a cărții Capitolul I Introducere, Funcții vectoriale și electrostatică O lecție, Natura, lasă-te să înveți despre tine Matthew Arnold Introducere În acest text subiectul calculului vectorial este prezentat în contextul electrostaticei simple Urmăm această procedură din două motive În primul rând, o mare parte din calculul vectorial a fost inventat pentru a fi utilizat în teoria electromagnetică și este ideal pentru acesta Prin urmare, această prezentare vă va arăta ce este calculul vectorial și, în același timp, vă va oferi o idee despre ce folosește În al doilea rând, avem convingerea profundă că matematica - în orice caz unele matematică - este cel mai bine discutată într-un context care nu este exclusiv matematic Astfel, vom pedala moale rigoarea matematică, Introducere, Funcții vectoriale, și electrostatică ceea ce credem că este un obstacol în calea învățării acestui subiect la o primă expunere la el și apelează pe cât posibil la intuiția fizică și geometrică Acum, dacă vrei să înveți calculul vectorial, dar știi puțin sau nimic despre electrostatică, nu trebuie să fii descurajat de abordarea noastră; nu sunt necesare cunoștințe foarte mari de fizică pentru a citi și înțelege acest text Sunt implicate doar cele mai simple caracteristici ale electrostaticei, iar acestea sunt prezentate în câteva pagini aproape de început Nu ar trebui să fie un impediment pentru nimeni De fapt, întreaga discuție se bazează pe căutarea unei metode convenabile de găsire a câmpului electrostatic având în vedere distribuția sarcinilor electrice care îl produc Acesta este firul care parcurge și unifică prezentarea noastră, astfel încât, cel puțin, tot ce trebuie să faceți este să ne credeți pe cuvânt pentru faptul că câmpul electric este o cantitate suficient de importantă pentru a justifica petrecerea timpului și a efortului în stabilirea unei metode generale de calcul al acesteia În acest proces, sperăm că veți învăța elementele de calcul vectorial După ce ați spus ceea ce nu trebuie să știți, acum trebuie să spunem ce trebuie să știți Pentru început, ar trebui, desigur, să fii fluent în calculul elementar În plus, ar trebui să știți cum să lucrați cu funcții ale mai multor variabile, derivate parțiale și integrale multiple (duble și triple) În cele din urmă, trebuie să știți ceva despre vectori Acesta este însă un subiect despre care prea mulți scriitori și profesori au făcut vreme grea Ceea ce ar trebui să știți despre el poate fi enumerat rapid: definiția vectorului, adunarea și scăderea vectorilor, multiplicarea vectorilor prin scalari, produse punctiforme și încrucișate și, în final, rezoluția vectorilor în componente O oră cu orice text rezonabil pe acest subiect ar trebui să vă ofere tot ce trebuie să știți despre el pentru a urma acest text Funcții vectoriale Una dintre cele mai importante cantități cu care ne ocupăm în studiul energiei electrice este câmpul electric și o mare parte din prezentarea noastră va folosi această cantitate Deoarece câmpul electric este un exemplu a ceea ce numim o funcție vectorială, începem discuția noastră cu un scurt rezumat al conceptului de funcție O funcție a unei variabile, scrisă în general у = f(x), este o regulă Ecuațiile diferențiale sunt utilizate într-o secțiune a acestui text Secțiunea nu este esențială și poate fi omisă dacă matematica este prea înspăimântătoare Funcții vectoriale care ne spune cum să asociem două numere x și y; dat x, funcția ne spune cum să determinăm valoarea asociată a lui y Astfel, de exemplu, dacă у = f(x) = x - , atunci calculăm у punând x la pătrat și apoi scăzând Deci, dacă x = , у = - = Funcțiile mai multor variabile sunt, de asemenea, reguli pentru asocierea unor seturi de numere De exemplu, o funcție de trei variabile desemnate w = F(x, y, z) spune cum să atribuiți o valoare lui w dat x, y și z Este util să privim acest concept geometric; luând (x, y, z) coordonatele carteziene ale unui punct din spațiu, funcția F(x, y, z) ne spune cum să asociem un număr cu fiecare punct Ca o ilustrare, o funcție T(x, y, z) poate da temperatura în orice punct (x, y, z) dintr-o cameră Funcțiile discutate până acum sunt funcții scalare; rezultatul "conectării" x în/(x) este scalarul у = f(x) Rezultatul "conectarii" celor trei numere x, y și z din T(x, y, z) este temperatura, un scalar Generalizarea la funcțiile vectoriale este simplă O funcție vectorială (în trei dimensiuni) este o regulă care ne spune cum să asociem un vector cu fiecare punct (x, y, z) Un exemplu este viteza unui fluid Desemnând această funcție v(x, y, z), ea specifică viteza fluidului, precum și direcția curgerii în punctul (x, y, z) În general, o funcție vectorială F(x, y, z) specifică o mărime și o direcție în fiecare punct (x, y, z) dintr-o regiune a spațiului Putem imagina o funcție vectorială ca o colecție de rânduri ar (Figura I- ), câte unul pentru fiecare punct (x, y, z) Direcția săgeții în orice punct este direcția specificată de funcția vectorială, iar lungimea acesteia este proporțională cu mărimea funcției Figura - O funcție vectorială, ca orice vector, poate fi rezolvată în compu- Introducere, Funcții vectoriale și electrostatică Figura - nente, ca în Figura - Lăsând i, j și к vectori unitari de-a lungul axelor x, y și, respectiv, z, scriem F(x, y, z) = iF/x, y, z) + jF/x, y, z) + kF,(x, y, z) Cele trei mărimi Fx, Fy și Fz, care sunt ele însele funcții scalare ale lui x, y și z, sunt cele trei componente carteziene ale funcției vectoriale F(x, y, z) într-un sistem de coordonate Un exemplu de funcție vectorială (în două dimensiuni pentru simplitate) este oferit de F(x, y) = ix + jy, care este ilustrat în Figura - Probabil recunoști asta Figura - Unii scriitori folosesc indicele pentru a indica derivata parțială; de exemplu, Fx = dF/dx Nu vom adopta o astfel de notație aici; indicele noștri vor desemna întotdeauna componenta vectorială Electrostatică funcţionează ca vector de poziţie r Fiecare săgeată din figură este în direcția radială (adică îndreptată de-a lungul unei linii care emană din fața originii) și are o lungime egală cu distanța sa față de origine Un al doilea exemplu, G(x, y) = ~iy + j* este prezentat în Figura - Ar trebui să verificați singuri că pentru această funcție vectorială toate săgețile sunt în direcția tangențială Figura - (adică fiecare este tangent la un cerc centrat la origine) și toate au aceeași lungime Electrostatică Ne vom baza discuția despre electrostatică pe trei fapte experimentale Primul dintre aceste fapte este existența sarcinii electrice în sine Există două tipuri de sarcină, pozitivă și negativă, și fiecare corp material conține sarcină electrică , deși adesea sarcinile pozitive și negative sunt prezente în cantități egale, astfel încât să existe o sarcină netă zero Al doilea fapt se numește legea lui Coulomb, după fizicianul francez care a descoperit-o Această lege prevede că electrostaticul Rețineți că prin convenție o săgeată este desenată cu coada, nu cu capul, în punctul în care este evaluată funcția vectorială Puriștii vor sublinia că neutronii, mezonii pi neutri, neutrinii și altele asemenea nu conțin încărcătură Introducere, Funcții vectoriale, și electrostatică ч Figura - forța dintre două particule încărcate (a) este proporțională cu produsul sarcinilor lor, (b) este invers proporțională cu pătratul distanței dintre ele și (c) acționează de-a lungul liniei care le unește Astfel, dacă qa și q sunt sarcinile a două particule aflate la o distanță r între ele (Figura - ), atunci forța pe q datorată lui q este F = K~-u, r unde u este un vector unitar (adică un vector de lungime ) care arată de la q la q , iar К este o constantă de proporționalitate În acest text vom folosi unități MKS raționalizate În acel sistem, lungimea, masa și timpul sunt măsurate în metri, kilograme și, respectiv, secunde, iar sarcina electrică în coulombs Cu această alegere a unităților К = /( тте ), unde constanta e , numită permisivitatea spațiului liber, are valoarea , X IO- coulombs per newton-metri , iar legea lui Coulomb se citește F = lună Q tt€ r ( - ) Ar trebui să vă convingeți că regula familiară "ca sarcinile se resping, spre deosebire de taxele se atrag" este inclusă în această formulă Al treilea și ultimul fapt se numește principiul suprapunerii Dacă F] este forța exercitată asupra qQ de qx când nu există alte sarcini în apropiere, iar F este forța exercitată asupra qt! de q când nu există alte sarcini în apropiere, atunci principiul suprapunerii spune că forța netă exercitată asupra qQ de q, și q când sunt ambele prezenți este suma vectorială Fj + F Aceasta este o afirmație mai profundă decât pare la prima vedere Se spune nu doar că forțele electrostatice se adaugă vectorial (toate forțele se adaugă vectorial), ci că forța dintre două particule încărcate nu este modificată de prezența altor particule încărcate Introducem acum o funcție vectorială a poziției, care va juca un rol principal în discuția noastră Este câmpul electrostatic, Electrostatică notat E(r) și definit prin ecuația E(r) = F(r)/g , sau F(r) = sunt constante (a) Cât de departe este obiectul de origine în orice moment t? (b) Găsiți viteza și accelerația obiectului în funcție de timp (c) Arătați că obiectul se mișcă pe calea eliptică - O sarcină + este situată în punctul ( , , ) și o sarcină - este situată în punctul (- , , ) Aflați câmpul electric al acestor două sarcini într-un punct arbitrar ( , y, ) pe axa y - În loc să folosim săgeți pentru a reprezenta funcții vectoriale (ca în problemele - și - ), uneori folosim familii de curbe numite linii de câmp O curbă у = y(x) este o linie de câmp a funcției vectoriale F(x, y) dacă în fiecare punct (x , y ) al curbei, F(x , y ) este tangentă la curbă (vezi figura ) Introducere, Funcții vectoriale și electrostatică (a) Să se arate că liniile de câmp у = y(x) ale unei funcții vectoriale F(x, y) = iFx(x, y) + jFv(x, y) sunt soluții ale ecuației diferențiale dy Fy(x, y) dx Fx(x,y)' (b) Determinați liniile de câmp ale fiecăreia dintre funcțiile problemei - Desenați liniile câmpului și comparați cu diagramele cu săgeți din problema - Capitolul II Integrale de suprafață și divergența O, aș putea să curg ca tine și să-ți fac pârâul Marele meu exemplu Sir John Denham Legea lui Gauss Deoarece câmpul electrostatic este o cantitate atât de importantă în electrostatică, rezultă că avem nevoie de o modalitate convenabilă de a-l găsi, având în vedere un set de sarcini La prima vedere s-ar putea părea că am rezolvat această problemă chiar înainte de a o afirma, căci, până la urmă, ecuațiile ( - ) și ( - ) nu ne oferă un mijloc de a găsi E? Răspunsul este, în general, nu Cu excepția cazului în care există foarte puține încărcări în sistem și/sau sunt aranjate simplu sau foarte simetric, suma din ecuația ( - ) și integrala din ecuația ( - ) sunt de obicei prohibitiv de dificile - și adesea impuse capabil - a efectua Astfel, aceste două ecuații oferă ceea ce este Integrale de suprafață si Divergenţă de obicei, doar o soluție "formală" a problemei, nu una practică și trebuie să cautăm o altă modalitate de a calcula câmpul E În cursul acestei turnări, ajungem inevitabil la acea relație remarcabilă cunoscută sub numele de legea lui Gauss Spunem "inevitabil" pentru că este greu să ne gândim la orice altă expresie în electricitate și magnetism elementar care conține câmpul electric [în afară, desigur, de ecuațiile ( - ) și ( - ), pe care le-am respins deja] Legea lui Gauss este І,Е'Й (II- ) Dacă vrem să încrucișăm toate t-urile și să punctăm toate z-urile, această integrală, strict vorbind, ar trebui scrisă ,y,z)' n(x, y, z) dS întrucât atât F cât și n sunt, în general, funcții de poziție Preferăm și, acolo unde este posibil, vom folosi notația mai puțin aglomerată F'iidS cu argumentele funcţiilor înţelese Suprafața S peste care integrăm o integrală de suprafață poate fi de două feluri: închisă sau deschisă O suprafață închisă, cum ar fi o înveliș sferică, împarte spațiul în două părți, una în interior și una în exterior, iar pentru a ajunge din interior în exterior, trebuie să treci prin suprafață O suprafață deschisă, cum ar fi o bucată de hârtie plată, nu are această proprietate; se poate ajunge dintr-o parte a foii pe cealalta fara a trece prin ea Definiția integralelor de suprafață dată în ecuația (II- ) se aplică la fel de bine atât suprafețelor închise, cât și celor deschise Cu toate acestea, integrala de suprafață nu este bine definită până când nu specificăm care dintre cele două direcții posibile ale normalei trebuie să folosim (vezi Figura II- ) În cazul unui Figura II- Afirmația "fiecare AS, " nu este destul de precisă Aria unui petic dreptunghiular, de exemplu ar putea tinde spre zero deoarece lățimea sa scade în timp ce lungimea rămâne fixă Acest lucru nu ar fi acceptabil Неге și în altă parte trebuie să interpretăm "fiecare AS, -> " pentru a însemna că toate dimensiunile liniare ale patch-ului /-lea tind să zero Integrale de suprafață si Divergenţă suprafață deschisă, direcția trebuie dată ca parte a enunțului problemei În cazul unei suprafețe închise, există un gentlemen's agreement care specifică direcția o dată pentru totdeauna: normalul este ales astfel încât să fie orientat spre exterior din volumul cuprins de suprafață Integrala din legea lui Gauss [Ecuația (II- )] este preluată pe o suprafață închisă Legea lui Gauss, de fapt, spune că integrala de suprafață a componentei normale a câmpului electric pe o suprafață închisă este egală cu sarcina totală (netă) cuprinsă de suprafață, împărțită la e Mai jos (paginile - și Problemele II- , - și - ) vom vedea cum, atunci când sarcinile sunt aranjate ordonat și simetric, legea lui Gauss poate fi folosită pentru a determina câmpul electric Dar scopul întregii noastre discuții va fi să supunem legea lui Gauss unei serii de aventuri chinuitoare care, în cele din urmă, o transformă într-o expresie utilă pentru a găsi E chiar și atunci când nu avem simetria care să ne ajute Uneori întâlnim integrale de suprafață care sunt puțin mai simple decât tipul pe care tocmai l-am definit, deși practic sunt aproape la fel Acestea sunt integrale de suprafață ale formei , y, z) dS, (П- ) unde integrând G(x, y, z) este o funcție scalară dată mai degrabă decât produsul scalar a două funcții vectoriale ca în (II- ) și (II- ) Noi definim acest tip de integrală de suprafață la fel cum am făcut mai sus: aproximăm S printr-un poliedru, formăm produsul G(xb yb zî) &Sb suma peste toate fețele și apoi luăm limita: N lim G(xz, yb Zi) AS, N->OO = fiecare AS -> (П- ) Ca exemplu de acest tip de integrală de suprafață, să presupunem că avem o suprafață de grosime neglijabilă cu densitatea suprafeței (adică masa pe unitate de suprafață) obținem masa totală a suprafeței Un exemplu de integrală de suprafață și mai simplă de acest fel este S Această integrală este luată ca definiție a suprafeței lui Evaluarea integralelor de suprafață Acum că am definit integralele de suprafață, trebuie să dezvoltăm metode pentru a le evalua și aceasta va fi sarcina noastră aici Pentru simplitate ne vom ocupa de integrale de suprafață de forma (П- ), unde integrând este o funcție scalară dată, mai degrabă decât forma puțin mai complicată (II- ) Nu va exista nicio pierdere a generalității în a face acest lucru, deoarece toate rezultatele noastre pot fi aplicate integralelor de forma (П- ) doar prin înlocuirea G(x, y, z) peste tot cu F(x, y, z) • n Pentru a evalua integrala s pe o porțiune S a suprafeței z = f(x, y) (vezi Figura II- ), ne întoarcem la definiția integralei de suprafață [Ecuația (II- )] Figura II- Integrale de suprafață si Divergenţă Figura - Strategia noastră va fi să reitim Д , zonei Д ?, a proiecției sale pe planul xy, așa cum se arată în Figura - Făcând acest lucru, așa cum vom vedea, ne va permite să exprimăm integrala de suprafață peste S în termenii unei integrale duble obișnuite peste R, care este proiecția lui S pe planul xy, așa cum se arată în Figura II- Relaționând A Z cu Д ?; nu este dificil dacă ne amintim că AS, (precum aria oricărei suprafețe plane) poate fi aproximată la orice grad de precizie dorit printr-un set de dreptunghiuri așa cum se arată în Figura - Figura - Din acest motiv trebuie doar să găsim relația dintre aria unui dreptunghi și proiecția lui pe planul xy Astfel, considerați un dreptunghi orientat astfel încât o pereche de laturi să fie paralelă cu planul xy (Figura - ) Dacă caii lungimile acestor laturi a, este clar că proiecțiile lor pe planul xy au și lungimea a Dar cealaltă pereche de laturi, de lungime b, au proiecții de lungime b' și, în general, b și b' nu sunt egale Astfel, pentru a reia aria dreptunghiului ab cu aria proiecției sale ab', trebuie să exprimăm Evaluarea Integrale de suprafață Figura - b în termeni de b' Acest lucru este ușor de făcut, deoarece dacă Ѳ este unghiul prezentat în figura П- , avem b = b'/cos Ѳ și astfel ab = ab' cos Ѳ ' Dacă lăsăm ii să desemnăm vectorul unitar normal dreptunghiului nostru, atunci ne putem convinge cu ușurință că cos Ѳ = n • к unde k, ca întotdeauna, este vectorul unitar în direcția z pozitivă Prin urmare, , ab' ab = -г n • к Deoarece aria ASZ poate fi aproximată cu o precizie arbitrară prin astfel de dreptunghiuri, rezultă că i\R, AS, = л n, • к unde, desigur, nz este vectorul unitar normal la suprafața plană Z Acum putem rescrie definiția integralei de suprafață [Ecuația (II- )] ca ff Л AK, G(x, y, z) dS = hm Уь Z/) л-~ JJs = n,*k fiecare AR -* (П- ) unde mențiunea "fiecare AS, -> " a fost înlocuită cu echivalent, dar acum mai potrivit "fiecare A ?z " Noi suntem Integrale de suprafață si Divergenţă acum, evident, pe drumul spre rescrierea integralei de suprafață peste S ca o integrală dublă peste R De fapt, lim ff n-"~ = П/'к J Jr n(x, у, z) *k fiecare A/?z-> (П- ) unde n(x, y, z) este vectorul unitar normal la suprafața S în punctul (x, y, z) Aceasta este o integrală dublă peste R, chiar dacă nu prea arată ca una Ceea ce pare să strice este acel z urât în G și n; o integrală dublă peste o regiune din planul xy în mod clar nu are nicio afacere care să conțină niciun z Dar dependența de z este falsă deoarece (x, y, z) sunt coordonatele unui punct pe S, și deci z = /(x, у) Prin urmare, în detrimentul de a face integrala să pară și mai acerbă decât este deja în ecuația (II- ), putem elimina dependența aparentă z a integrării și scriem л '-' -' dx dy rn[x, y,/(x, y)] • к (П- ) Cei slabi de inimă pot lua curaj; în cele mai multe cazuri, această integrare se reduce rapid la ceva mult mai simplu și mai plăcut, fapt pe care îl vom demonstra prin exemplu mai jos În acest moment introducem expresia pentru vectorul normal unitar [Ecuația (II- )] Găsim Vi + (а//Эх) + (Э/ Эу) ' și astfel ecuația (II- ) devine y, z) dS dxdy (П- ) Astfel, integrala de suprafață a lui G(x, y, z) peste suprafața S a fost exprimată ca o integrală dublă a unei funcții cu aspect dezordonat peste regiunea R, proiecția lui în planul xy După cum am remarcat mai sus, în practică integrala este de obicei mult mai puțin îngrozitoare decât ea Evaluarea Integrale de suprafață apare scris în ecuațiile (II- ) sau ( - ) Veți vedea acest lucru în exemplele pe care le oferim acum Să calculăm mai întâi integrala de suprafață + z) dS unde S este porțiunea planului x + у + z = în primul oc-tant prezentat în Figura II- (a) Proiecția lui pe planul xy este triunghiul R prezentat în figură Ecuația lui S poate fi scrisă z = /(*,?) = - x - у din care obținem У dx dy astfel încât Prin urmare + z)dS = (x + z) dx dy = + - x - y) dxdy ( - y) dx dy, Integrale de suprafață si Divergenţă unde am folosit z = - x - y Aceasta este o integrală dublă simplă cu valoarea l/Ѵз, așa cum ar trebui să puteți verifica Ca un al doilea exemplu, să calculăm integrala de suprafață unde S este octantul sferei cu raza centrată la origine, așa cum se arată în Figura II- (b) Proiecția lui pe Figura П- (b) planul xy (adică R) este aria cuprinsă de sfert de cerc Ecuația lui S este x + y + z = , sau z =Жу) = +V - x - у Rezultă atunci că Эх z și dy Z' astfel încât = lv + +?-l unde am folosit x + y + z = Prin urmare, z dx dy - jj zdxdy Evaluarea Integrale de suprafață Înlocuind z în termeni de x și y, obținem /^' Aceasta este o integrală dublă obișnuită și ar trebui să verificați dacă valoarea sa este rr/ [Sugestie: Convertiți în coordonate polare: x - r cos Ѳ și у = r sin Ѳ Integrarea este atunci banală ] Trebuie subliniat că discuția de mai sus sa bazat pe presupunerea că suprafața S este descrisă printr-o ecuație de forma z = f(x, y); într-o astfel de situație o integrală de suprafață este convertită într-o integrală dublă peste o regiune din planul xy Dar se poate întâmpla ca o suprafață dată să fie descrisă mai convenabil printr-o ecuație de forma у = g(x, z) ca în Figura II- (a) Dacă este așa, atunci y, z) dS Figura II- (a) unde R este o regiune în planul xz În mod similar, dacă avem o suprafață descrisă de x = h(y, z), ca în Figura II- (b), atunci folosim s , y, z) dS eu dy dz, R Integrale de suprafață si Divergenţă Figura П- (Ь) unde R în acest caz este o regiune în planul yz În cele din urmă, o suprafață poate avea mai multe părți și apoi poate fi convenabil să proiectați diferite părți pe planuri de coordonate diferite Pentru a evalua integralele de suprafață de forma (II- ), adică pur și simplu înlocuim G cu F • n în ecuația ( - ) pentru a obține /// • dS - Я"F • v + (I) + (Ю * Dacă folosim acum ecuația (II- ) pentru a scrie acest lucru în detaliu, aflăm că factorul rădăcină pătrată se anulează și obținem - ^ [x, у, f(x, у)] + F[x, y,f(x, у)] l dx dy (П- ) Lăsăm cititorului să scrie formulele analoage atunci când suprafața este dată de у = g(x, z) sau x = h(y, z), care trebuie proiectat pe regiuni din xz- și yz- avioane, respectiv Această ultimă ecuație [Ecuația ( - )] este suficientă pentru a face bărbații puternici să plângă, dar, ca și înainte, în majoritatea calculelor se reduce rapid la ceva destul de blând De exemplu, să presupunem că dorim să calculăm ffsF • n dS, unde F(x,y, z) = iz ~ jy + kx și Sis proporția planului x + y + z = Evaluarea Integrale de suprafață z Figura II- (a) delimitat de planurile de coordonate, adică triunghiul care se înclină grațios în Figura II- (a) Vectorul normal n este ales astfel încât să fie orientat departe de origine, așa cum se arată în Figura II- (a), și vom proiecta pe planul xy Avem Z=/(x,y)= -jy, și așa dx ' dy De asemenea avem = z = - | - y, Fy= -y, F =x Prin urmare (~l) + y(~i) + x dx dy > , ] , Regiunea R peste care trebuie luată integrala este prezentată în Figura II- (b) Problema a fost astfel redusă la сотри- Integrale de suprafață si Divergenţă Figura II- (b) integrarea unei integrale duble destul de simple și ar trebui să realizați singur integrarea (răspunsul este |) Ca un al doilea exemplu, să presupunem F(x, y, z) =ixz + kz Pentru S să luăm octantul suprafeței sferice prezentate în figură - Atunci z = f(x,y) = Vi - x~ - y, anj am arătat deja (vezi pagina ) că = x dx ? și , volumul inclus Д V trebuie, desigur, să se apropie și el Divergenta Figura - zero Astfel, fluxul tinde, de asemenea, spre zero și afirmația este dovedită Nu numai că am dat o dovadă, dar (și acesta este ideea) acum putem izola o mărime care nu dispare ca S Ecuația de împărțire ( - ) prin Д V, obținem E • n dS - Рду/^д Această expresie, oricât de incomodă și neatrăgătoare ar fi, este totuși apropiată de ceea ce căutăm, deși implică încă o integrală a lui E pe o întreagă suprafață Căci dacă luăm acum limita ca se micșorează la zero în jurul unui punct din Д V ale cărui coordonate sunt (x, y, z), atunci, după cum vedem din ecuația ( - ), densitatea medie рдѵ se apropie de p(x , y, z), densitatea la (x, y, z) și obținem lim лѵ f IE ' " dS = p(x'y' z^/e°- (П- ) дѵ->о ілѴ J Js despre (x,y,z) Această expresie este, desigur, de-a dreptul hidosă, și dacă va fi de vreun folos practic, depinde de faptul că suntem capabili să punem partea stângă într-o formă care arată și acționează cel puțin pe jumătate civilizat Ne întoarcem la această sarcină acum Divergenta Să considerăm integrala de suprafață a unei funcții vectoriale arbitrare F(x, y, z): Această linie de raționament și concluzia trebuie modificate dacă sistemul conține taxe punctuale Integrale de suprafață si Divergenţă Ne va interesa raportul dintre această integrală și volumul cuprins de suprafața S pe măsură ce volumul se micșorează la zero aproximativ un punct, pentru că acesta este exact tipul de mărime care apare în ecuația ( - ) Această limită este suficient de importantă pentru a garanta un nume și o notație specială Se numește divergența ofF și este denumită div F Astfel, div F = lim -J- Г f F • n dS ( - ) дѵ-"o AVJ Js despre (x,y,z) Această cantitate este în mod clar un scalar Mai mult, va avea, în general, valori diferite în puncte diferite (x, y, z) Astfel, divergența unei funcții vectoriale este o funcție scalară de poziție Ecuația ( - ) poate fi acum scrisă div E = p/e ( - ) În această etapă, totuși, noua noastră notație fantezică are doar o valoare cosmetică, ajutând la înfrumusețarea unei ecuații urâte Dacă are și vreo valoare practică este problema abordată în următoarea discuție, în care calculăm de fapt limita raportului fluxului la volumul închis și aflăm că poate fi exprimat în mod rezonabil simplu în termeni de anumite derivate parțiale Înainte de a trece la acest calcul, totuși, merită menționat că, dacă luăm noua noastră terminologie la propriu, putem interpreta ecuația ( - ) în sensul că câmpul "diverge" de la un punct și cât de mult diverge, ca să spunem așa , depinde de cât de multă sarcină există în acel punct, așa cum este reprezentată de densitatea de acolo Următoarea noastră ordine de lucru este să găsim expresia relativ simplă pentru divergența unei funcții vectoriale promisă mai sus Astfel, luați în considerare un mic paralelipiped dreptunghiular cu muchii de lungime Ax, Ay și Az paralele cu axele de coordonate (Figura - ) Fie punctul din centrul micului cuboid să aibă coordonatele (x, y, z) - Calculăm integrala de suprafață a lui F peste suprafața cuboidului considerând integrala ca o sumă de șase termeni, câte unul pentru fiecare cuboid față Începem prin a lua în considerare fața marcată A] din figură Noi vrem De acum înainte ne vom referi la acesta drept "cuboid", un termen inventat care necesită mai puțin timp și spațiu decât "paralepipedul dreptunghiular" sesquipedal Divergenta Figura - Acum este clar că vectorul unitar normal la această față și îndreptat spre exterior din volumul inclus este i Astfel, deoarece F • i = Fx, integrala precedentă este y, z) dS Prin presupunere, cuboidul este mic (în cele din urmă vom lua limita pe măsură ce se micșorează la zero) Prin urmare, putem calcula această integrală aproximativ ca Fx evaluată la centrul feței Sj înmulțit cu aria feței Coordonatele centrului lui Sx sunt (x + Дх/ , y, z) Prin urmare, FJx + ^,y,zl ДуДг (П- ) Raționalul din spatele acestui lucru este următorul: Există o teoremă a valorii medii, care ne spune că integrala lui F, peste , este egală cu aria lui , înmulțită cu funcția evaluată undeva pe S, Deoarece S, este mic, punctul în care ar trebui să evaluăm Fx și punctul în care îl evaluăm (adică centrul) trebuie să fie aproape unul de celălalt, iar Fx trebuie să aibă aproape aceeași valoare în cele două puncte Prin urmare, ceea ce ne oferă procedura noastră este o bună aproximare a valorii integralei Mai mult, pe măsură ce cuboidul se micșorează la zero, cele două puncte se apropie din ce în ce mai mult, astfel încât în limită rezultatul nostru [Ecuația (П- )] va fi exact Integrale de suprafață si Divergenţă Același tip de raționament aplicat feței opuse S [a cărei normală exterioară este - i și al cărei centru este la (x - Дх/ , y, z)] conduce la = - Fj x - > У> z I Ду Az (П- ) Adunând contribuțiile acestor două fețe [Ecuațiile ( - ) și ( - )], obținem Recunoscând că Ax Ду Дг = Д V, volumul cuboidului, avem ff F • n dS ДѴ J Js,+s ( - ) Acum trebuie să luăm limita acesteia pe măsură ce ДѴ se apropie de zero Dar, desigur, pe măsură ce ДѴ ajunge la zero, la fel face fiecare dintre laturile cuboidului Astfel, în partea dreaptă a ecuației ( - ) putem scrie lim^^o în locul limAV " și găsim lim ~ jf F • n dS av-o ДѴ J JS +S, р/ Дх [ " I Ax ] F^x + T,y,z^-F^xT,y,zj dF = llm -л = Дх-" Дх ох Rețineți că am amânat calcularea contribuțiilor de la celelalte patru fețe ale cuboidului Divergența evaluată la (x, y, z) Această ultimă egalitate rezultă din definiția derivatei parțiale Nu ar trebui să fie surprinzător faptul că celelalte două perechi de fețe ale cuboidului contribuie cu dFy!dy și dF Jdz Prin urmare, ff л dFx dF dF lim - F • n dS = - + - - + -r-Л ЛѴ-" Д VJ Js ax dy dz Limita din partea stângă a acestei ultime ecuații, așa cum am remarcat deja, este divergența lui F [Ecuația ( - )] Astfel, tocmai am demonstrat asta dF dFv dF di v F = -^ + - + ( - ) dx dy dz Se poate arăta că acest rezultat este independent de forma volumului folosit pentru a-l obține (vezi problema - ) Utilizarea ecuației ( - ) pentru a găsi divergența unei funcții vectoriale este o chestiune simplă, dar vom da un exemplu doar pentru înregistrare Luați în considerare funcția F(x, y,z) = ix + jxy + kyz Avem dx x' dy dF și -= y dz J Prin urmare, div F = x + x + у = x + y Revenind acum la câmpul electrostatic, combinăm ecuațiile ( - ) și ( - ) pentru a obține dEx dey dE аГ + IF + FF ~ p/e°- ( - ) Această ecuație, care este mult mai generală decât sugerează derivarea noastră, este una dintre ecuațiile lui Maxwell și este complet echivalentă cu legea lui Gauss [Ecuația (П-l)] Este uneori numită "forma diferențială" a legii lui Gauss Am ajuns acum la scopul nostru (aproape!), deoarece am legat o proprietate a câmpului electrostatic într-un punct (adică divergența acestuia) cu o cantitate cunoscută (densitatea de sarcină) în acel punct Cu toate acestea Integrale de suprafață si Divergenţă trebuie spus că Ecuația ( - ) poate fi privită într-un sens ca o singură ecuație (diferențială) în trei necunoscute (E" Ey, E,) și din acest motiv nu este adesea folosită în această formă pentru a găsi câmpul Totuși, rezultă că cele trei componente ale lui E pot fi legate între ele foarte elegant; când vom dezvolta această relație, ne vom întoarce la această întrebare de a găsi un mijloc convenabil de a calcula E Divergența în coordonatele cilindrice și sferice Se vede adesea ecuația ( - ) dată ca definiție a divergenței funcției vectoriale F Deși acest lucru este cu siguranță acceptabil, preferăm mult să definim divergența ca limită a fluxului la volum așa cum este menționat în ecuația ( ) - ) Ecuația ( - ) este atunci doar forma pe care o ia divergența în coordonate carteziene În alte sisteme de coordonate arată destul de diferit De exemplu, în coordonatele cilindrice, funcția F are trei componente, pe care nu veți fi șocați să le învățați sunt desemnate Fr, Fe și F [vezi Figura II- (a)] Pentru a obține divergența lui F în cilindru Figura П- (а) coordonatele drice, considerăm "cuboidul cilindric" prezentat în Figura II- (b) cu volum AV = r Ar ДѲ Az și centru în punctul (г, Ѳ, z) Fluxul lui F prin fața marcată este F • n dS = - Fr I r + Ѳ, z II r + I ДѲ Az, Rețineți că în cazul cartezian (Figura - ) fiecare față a cuboidului este dată de o ecuație de forma x = constantă, у = constantă sau z = constantă În același mod, fiecare față a suprafeței din Figura II- (b) este dată de o ecuație de forma r = constantă, Ѳ = constantă sau z = constantă Divergența în coordonatele cilindrice și sferice Figura D- (b) în timp ce prin faţa marcată este Т'ѳ' )(''~т)АвЛг- Adunând aceste două rezultate și împărțind la volumul А У al cuboidului, găsim n dS r Ar Ar Ar г I Ar Q Frlr-T, ,Z care în limita pe măsură ce Ar (şi deci AV) se apropie de zero devine Argumentând în mod analog pentru celelalte patru fețe (vezi problema - ), ajungem în sfârșit la expresia divergenței în coordonate cilindrice: până la dFe div F (rF) - - rarг эѳ az ( - ) , rx râ?W Integrale de suprafață si Divergenţă Figura П- În coordonatele sferice în care componentele lui F sunt Fr, F și F$ (vezi Figura П- ), un raționament similar (vezi Problema - ) conduce la expresia div F = - ~ (HFr) H ; - (sin фГф) ; т - r дг г sin ф Эф ф г sin ф ЭѲ (П- ) Notația Del Există o notație specială în care poate fi scrisă divergența Ar fi puțin sau deloc motiv pentru introducerea acestuia dacă ar servi doar pentru a oferi un alt mod de a scrie "div", dar, după cum vom vedea în curând, are o utilitate considerabilă în calculul vectorial Să deținem o cantitate desemnată V (a se citi "del") prin următoarea ecuație destul de ciudată: w д , Э , Э V = - + j v- + k - dx dy dz Dacă luăm produsul scalar al lui V și o funcție vectorială F = iFx + jFv + kF , obținem V"F= i|- + j|- + k/- • \ dx dy dz) = ^-F + Ip +-F dx x dy y dz z (iFx + jFy + кЯ) Acum interpretăm "produsul" lui д/dx și Fx ca o derivată parțială; acesta este, Э dF, - F = dx x dx Teorema divergenței Există ecuații similare pentru celelalte două "produse" (d/dy)Fy și (d/dz)Fv Cu această convenție recunoaștem V • F ("del dot F") ca fiind la fel ca div F și, de acum înainte, la conform practicii moderne de notație, vom folosi întotdeauna V • F pentru a indica divergența Astfel, se vor scrie ecuațiile ( - ) și ( - ) V • E = p/e Matematicienii caii un simbol ca V un operator Când "operăm" cu V punând-o într-o funcție vectorială, obținem divergența acelei funcție, așa cum tocmai am văzut În discuțiile ulterioare vom introduce alte trei mărimi (gradient, curl și laplacian), toate fiind operatori și toate putând fi scrise în termeni de V Teorema divergenței Pentru restul acestui capitol, ne abatem de la curentul principal al narațiunii noastre pentru a discuta o teoremă faimoasă care afirmă o legătură remarcabilă între integralele de suprafață și integralele de volum Deși această relație poate fi sugerată de munca pe care am făcut-o în electrostatică, teorema este o afirmație matematică care se susține în circumstanțe destul de generale Este independent de orice fizică și este aplicabil în multe locuri diferite Se numește teorema divergenței și uneori teorema lui Gauss (a nu se confunda cu legea lui Gauss) Nu vom oferi o demonstrație riguroasă din punct de vedere matematic a teoremei divergenței; o asemenea dovadă este dată în multe texte de calcul avansat În schimb, prezentăm aici o dovadă brută și gata a unui alt fizician Astfel, luați în considerare o suprafață închisă S Subdivizați volumul V închis de S în mod arbitrar în N subvolume, dintre care unul este prezentat în Figura - (desenat ca un cub pentru comoditate) Noi Figura - Integrale de suprafață si Divergenţă începeți demonstrația noastră afirmând că fluxul unei funcții vectoriale arbitrare F(x, y, z) prin suprafața S este egal cu suma fluxurilor prin suprafețele fiecăruia dintre subvolume: E ff F"ndS (= J -Is, ( - ) Неге S[ este suprafața care cuprinde subvolumul Д V, Pentru a stabili ecuația ( - ), luați în considerare două subvolume adiacente (Figura - ) Fie ca chipul lor comun să fie notat Deci Fluxul prin Figura - subvolumul marcat cu în Figura П- include, desigur, o contribuție din So, care este ăj dS Неге n, este un vector unitar normal cu fața So și, prin convenția noastră obișnuită, indică spre exterior din subvolumul Fluxul prin subvolumul marcat include și o contribuție din So: | F"n Punând împreună ecuațiile de la ( - ) la ( - ), ajungem la rezultatul nostru: V • F dV ( - ) Aceasta este teorema divergenței Cu cuvinte, se spune că fluxul unei funcții vectoriale printr-o suprafață închisă este egal cu triplul inte- Integrale de suprafață si Divergenţă grai al divergenței acelei funcții asupra volumului închis de suprafață Motivul principal pentru care demonstrația dată mai sus nu este riguroasă este că o integrală triplă este definită ca limita unei sume de formă S(xb yb Zi) A Vb i unde funcţia g este bine definită În ecuația ( - ), totuși, cantitatea care înmulțește elementul de volum AV; în fiecare termen al sumei nu este o funcţie bine definită în acest sens Adică, deoarece AVZ tinde să zero, cantitatea dintre paranteze pătrate se modifică; poate fi identificat ca divergenta lui F numai in limita Un tratament atent și riguros ar arăta că ecuația ( - ) este validă dacă F (adică Fx, Fy și F ) este continuă și diferențiabilă, iar derivatele sale primare sunt continue în V și pe S Acum să ilustrăm teorema divergenței Deoarece pagini nesfârșite de integrale hidoase nu ne vor servi scopului, vom folosi un exemplu simplu Lăsa F(x, y, z) = ir + jy + kz și alegeți pentru suprafața prezentată în Figura - , constând din învelișul emisferic cu raza și regiunea R a Figura - planul xy închis de cercul unitar Pe emisferă avem n = ix + jy + kz, astfel încât n*F=x + y + z = Astfel, pe emisferă, JJF • ndS = JJ dS = ir, Două aplicații simple ale teoremei divergenței unde ultima egalitate rezultă din faptul că integrala este doar aria suprafeței emisferei unității Pe regiunea R avem ii = -к astfel încât n • F = -z Prin urmare, pe R, F • n dS = deoarece z = peste tot pe R Astfel, nu există nicio contribuție la integrala de suprafață din regiunea circulară R și DACĂ • n dS = ir În continuare aflăm printr-un calcul banal că V • F = Rezultă atunci că unde folosim faptul că volumul emisferei unitare este -тт/З Deoarece integralele de suprafață și volum sunt egale, aceasta ilustrează ecuația ( - ) Două aplicații simple ale teoremei divergenței Ca exemplu de utilizare a teoremei divergenței, oferim o derivație alternativă a ecuației (П- ), a cărei analiză ne-a condus la teorema divergenței în sine Cu alte cuvinte, atât de ușor ar fi fost dacă am fi cunoscut teorema divergenței de la început! Începem cu legea lui Gauss în formă Ле,й" ЧЛЛл ( - ) Apoi aplicăm teorema divergenței integralei de suprafață din ecuația de mai sus pentru a obține E • ii dS = v ( - ) Astfel, combinând ecuațiile (П- ) și ( - ), găsim Integrale de suprafață si Divergenţă În general, dacă două integrale de volum sunt egale, nu este neapărat adevărat că integralele lor sunt egale S-ar putea ca integralele să fie egale numai pe volumul particular al integrării V, iar prin integrarea pe un volum diferit, am distruge egalitatea În cazul de față, totuși, acest lucru nu este adevărat deoarece legea lui Gauss este valabilă pentru orice volum arbitrar V și nu putem deranja egalitatea schimbând volumul Dar acest lucru poate fi așa numai dacă integranții sunt egali Prin urmare, V • E = p/e , care ar trebui să pară cunoscut! Un alt exemplu de utilizare a teoremei divergenței este următorul Să presupunem că într-o regiune a spațiului "lucrurile" (materie, sarcină electrică, orice) se mișcă (Figura - ) Lasă bârlogul- Figura - Siitatea acestor lucruri în orice punct (x, y, z) și în orice moment t să fie p(x, y, z, f) și să fie viteza v(x, y, z, t) Mai mult, să presupunem că aceste lucruri sunt conservate; adică nu este nici creată, nici distrusă Concentrându-ne pe un volum arbitrar V din spațiu, ne întrebăm: Care este viteza cu care se schimbă cantitatea de lucruri din acest volum? În orice moment t cantitatea de lucruri din V este jj jp(-*'y> z> dv iar ritmul cu care se schimbă este Două aplicații simple ale teoremei divergenței (Pentru a putea muta derivata sub semnul integral în acest fel, este nevoie ca dp/dt să fie continuu ) În continuare ne amintim dintr-o discuție anterioară că viteza cu care materialul curge printr-o suprafață S este | pv • n dS Apoi afirmăm că viteza cu care se schimbă cantitatea de material din V este egală cu viteza cu care curge prin suprafața de închidere ; sub formă de ecuație se citește această afirmație V | pv • n dS Există două caracteristici despre această ecuație care necesită discuție: Semnul negativ trebuie inclus deoarece integrala de suprafață, așa cum este definită, este pozitivă pentru un debit net din volum, dar un flux net de ieșire înseamnă că cantitatea de material din volum este în scădere Această ecuație afirmă că cantitatea de obiecte din V se poate modifica numai ca urmare a fluxului de lucruri peste granița Dacă lucrurile ar fi create sau distruse în V, termenii ar trebui să fie incluși în ecuație pentru a reflecta acest fapt Absența oricăror astfel de termeni este astfel o expresie a conservării materialelor Acum, în sfârșit, să aplicăm teorema divergenței Găsim | pv • n dS = Prin urmare, v Argumentând, așa cum am făcut mai sus, că V este un volum arbitrar, putem spune atunci ^=-V-(Pv) ( - ) Integrale de suprafață De obicei definim densitatea de curent J = pv și scriem ecuația și ( - ) ca Divergenţă Ș + VJO O ecuație de acest tip este denumită ecuație de continuitate și este, după cum am văzut, o expresie a unei legi de conservare (vezi Problemele III- , Ш- și IV- ) Pe lângă faptul că joacă un rol important în teoria electromagnetică, este o ecuație de bază atât în hidrodinamică, cât și în teoria difuziei În sfârșit, în analiza fluxului de căldură sunt implicate considerații similare cu cele care au condus la ecuația de continuitate PROBLEME II-l Găsiți un vector unitar n normal pentru fiecare dintre următoarele suprafețe (a) z = - x - y (c) z = ( - x ) / (b) z = (x + y ) / (d) z = x + y (e) z = ( - x /a - y /a )' - (a) Să se arate că vectorul unitar normal pe plan ax + by + cz = d este dat de n = ±(ia + jb + kc)/(a + b + c ) / (b) Explicați în termeni geometrici de ce această expresie pentru n este independentă de constanta d II- Deduceți expresii pentru vectorul normal unitar pentru suprafețe date de у = g(x, z) și de x = h(y, z) Folosiți fiecare pentru a rededuce expresia normalei la planul dat în problema II- II- În fiecare dintre următoarele, utilizați ecuația - pentru a evalua integrala de suprafață Jf , G(x, y, z) dS (a) G(x, y, z) = z, unde este porțiunea planului x + у + z = în primul octant (b) G(x,y,z)= j- + (x + y ) unde este porțiunea paraboloidului z = x + y între z = și z = (c) G(x,y,z) = (l-x -y )* unde S este emisfera z = ( - x - y ) / Probleme II- În fiecare dintre următoarele, utilizați ecuația П- pentru a evalua integrala de suprafață f fs F • n dS (a) F(x, y, z) = ix - kz, unde S este porțiunea planului x + у + z = în primul octant (b) F(x, y, z) = ix + jy + kz, unde este emisfera z = Va - x - y (c) F(x, y, z)=yy + k, unde este porțiunea paraboloidului z = - x - y deasupra planului xy - Distribuția masei pe învelișul emisferic Z = (Л - x - y ) ' este dat de а(х,у,г) = (а /Я )(х + у ) unde cr este o constantă Găsiți o expresie în termeni de ct și R pentru masa totală a învelișului - Aflați momentul de inerție în jurul axei z a învelișului emisferic a problemei П- - Un câmp electrostatic este dat de E = X(iyz + jxz + kxy), unde X este o constantă Folosiți legea lui Gauss pentru a găsi sarcina totală înconjurată de suprafața prezentată în figură constând din ,, emisferă Z = (Л - x - y ) ' , și , baza sa circulară în planul xy - Un câmp electrostatic este dat de E = X(ix + jy), unde X este o constantă Folosiți legea lui Gauss pentru a găsi sarcina totală cuprinsă de suprafața prezentată în figură constând din S ; porţiunea curbată a semicilindrului z = (r - y )l/ de lungime h\ S şi , cele două piese de capăt plane semicirculare; și S , porțiunea dreptunghiulară a planului xy Exprimați rezultatele în termeni de X, r și h Integrale de suprafață si Divergenţă - Se întâmplă uneori ca integralele de suprafață să poată fi evaluate fără a utiliza procedurile de lungă durată prezentate în text Încercați să evaluați ffs F • n dS pentru fiecare dintre următoarele; gândește-te puțin și evită multă muncă! (a) F = ix + jy + kz S, cele trei pătrate fiecare cu latura b, așa cum se arată în figură (b) F = (ix + jy) ln(x + y ) S, cilindrul (inclusiv partea de sus și de jos) cu raza Л și înălțimea h prezentate în figură X Probleme (c) F = (ix + jy + kz)e '(?+yi+?) S, suprafața sferei cu raza R centrată la origine, așa cum se arată în figură (d) F = i£(x), unde E(x) este o funcție scalară arbitrară a lui x S, suprafața cubului cu latura b prezentată în figură - (a) Folosiți legea lui Gauss și simetria pentru a găsi câmpul electrostatic în funcție de poziție pentru un plan de sarcină uniform infinit Lăsați sarcina să se afle în planul yz și notați sarcina pe unitate de suprafață cu sau (b) Repetați partea (a) pentru o placă infinită de sarcină paralelă cu planul yz, a cărei densitate este dată de p(x) = Po , -b Z>, unde p și b sunt constante (c) Repetați partea (b) cu p(x) = poe ^/b' - (a) Folosiți legea lui Gauss și simetria pentru a găsi câmpul electrostatic în funcție de poziție pentru o linie de sarcină uniformă infinită Lasă sarcina să se afle de-a lungul axei z și notăm sarcina pe unitatea de lungime cu X Integrale de suprafață si Divergenţă (b) Repetați partea (a) pentru un cilindru infinit de sarcină a cărui axă coincide cu axa z și a cărui densitate este dată în coordonate cilindrice de ₽w - {*• r b unde p și b sunt constante (c) Repetați partea (b) cu p(r) = poe - (a) Folosiți legea lui Gauss și simetria pentru a găsi câmpul electrostatic în funcție de poziție pentru distribuția de sarcină simetrică sferic a cărei densitate este dată în coordonate sferice de pfr) = Po , r b, unde p și b sunt constante (b) Repetați partea (a) pentru p(r) = poe (c) Repetați partea (a) pentru P(r) = r b Cum trebuie să fie legate p și p, astfel încât câmpul să fie zero pentru r > b? Care este taxa totală a acestei distribuții în aceste circumstanțe? - Calculați divergența fiecăreia dintre următoarele funcții folosind ecuația ( - ): (a) ix + jy + kz (b) iyz + jxz + kxy (c) ie~x + je~y + ke~z (d) i - j + kz (e) (~ixy + jx )/^ + y ), (x, y) * ( , ) (f) kVx + y (g) ix + jy + kz (h) (-iy + jx)/Vx + y , (x, у) * ( , ) - (a) Calculați ffs F • n dS pentru funcția din problema II- (a) pe suprafața unui cub cu latura s al cărui centru este la (x , y , Zo) și ale cărui fețe sunt paralele cu planuri de coordonate (b) Împărțiți rezultatul de mai sus la volumul cubului și calculați limita coeficientului ca s -" Comparați rezultatul dvs cu divergența găsită în problema П- (а) (c) Repetați părțile (a) și (b) pentru funcția Problemei II- (b) și (c) - (a) Calculați divergența funcției F(x, y, z) = i/(x) + j/(y) + k/(- z) și arată că este zero în punctul (c, c, -c/ ) Probleme (b) Calculați divergența lui G(x, y, z) = dacă(y, z) + jg(x, z) + кй(х, у) II-l În text am obţinut rezultatul " " FZ dFy r)F V'F = ^ + - + ~H dx oy dz prin integrarea pe suprafaţa unui mic paralelipiped dreptunghiular Ca exemplu al faptului că acest rezultat este independent de suprafață, rederivați-l folosind suprafața în formă de pană prezentată în figură - (a) Fie i, j și к vectori unitari în coordonate carteziene și er, ее și e, vectori unitari în coordonate cilindrice Arata asta i = er cos Ѳ - ёѳ sin Ѳ, j = ё, sin Ѳ + ёѳ cos Ѳ , к = ёг (b) Rescrieți funcția din problema II- (e) în coordonate cilindrice și calculați divergența acesteia, folosind ecuația ( - ) Convertiți rezultatul înapoi în coordonate carteziene și comparați cu răspunsul obținut în problema П- (е) (c) Repetați partea (b) pentru funcția Problemei II- (f) - (a) Fie i, j și к vectori unitari în coordonate carteziene și er, ёѳ și ёф vectori unitari în coordonate sferice Arata asta i = ёг sin ф cos Ѳ + ёф cos ф cos Ѳ - ёѳ sin Ѳ j = ёг sin ф sin Ѳ + ёф cos ф sin Ѳ + ёѳ cos Ѳ к = ёг cos ф - ёф sin ф [Sugestie: Este mai ușor să exprimați ё" ёв și ёф în termeni de i, j și к și apoi rezolvați algebric pentru i, j și к Pentru a face acest lucru, utilizați mai întâi faptul că er - rlr = (ix + jy + kz)/r În continuare, raționând geometric, arătați că ёѳ = -i sin Ѳ + j cos Ѳ În cele din urmă, calculați ёф = ёѳ X er J (b) Rescrieți funcția problemei II- (g) în coordonate sferice și calculați divergența acesteia folosind ecuația (П- ) Convertit Integrale de suprafață si Divergenţă rezultatul înapoi la coordonatele carteziene și comparați cu răspunsul obținut în problema II- (g) (c) Repetați partea (b) pentru funcția Problemei II- (h) - În coordonatele cilindrice divergenţa lui F este dată de V F = - - (rF ) + - - + r drv r) r ЭѲ dz În text (paginile -ЛЗ) am derivat primul termen al acestei expresii Procedând la fel, obțineți ceilalți doi termeni П- Repetați problema - pentru a obține divergența în coordonate sferice realizând integrala de suprafață pe suprafața volumului prezentat în figură și obținând astfel expresia vF=M(r Fr)+ până la dFe (sinW + -? - Se consideră o funcție vectorială de forma F(r) = ej(r), unde er = (ix + jy + kz)/r este vectorul unitar în direcția radială, r = (x + y + z ) / și f(r) este o funcție scalară diferențiabilă Folosind rezultatele problemei - , determinați /(r) astfel încât V • F = Se spune că o funcție vectorială a cărei divergență este zero este solenoidală - Verificați teorema divergenței în fiecare dintre următoarele cazuri: (a) F = ix + jy + kz S, suprafața cubului cu latura b prezentată în figură Probleme (b) F = err + ez, r = ix + jy S, suprafața sfertului de cilindru (raza R, înălțimea h) prezentată în figură (c) F = er r , r = ix + jy + kz S, suprafața sferei cu raza R centrată la origine, așa cum se arată în figură Integrale de suprafață Ц- și Divergenţă (a) Una dintre ecuațiile lui Maxwell afirmă că V • В = , unde В este orice câmp magnetic Arata asta В • n dS = pentru orice suprafață închisă S (b) Determinați fluxul unui câmp magnetic uniform В prin suprafața curbată a unui con circular drept (raza R, înălțimea h) orientat astfel încât В să fie normal la baza conului, așa cum se arată în figură (Un câmp uniform este unul care are aceeași mărime și direcție peste tot ) - Folosiți teorema divergenței pentru a arăta că n dS = , unde S este o suprafață închisă și n este vectorul unitar normal la suprafața П- (a) Folosiți teorema divergenței pentru a arăta că unde S este o suprafață închisă care cuprinde o regiune de volum V, n este un vector unitar normal cu suprafața S și r = ix + jy + kz (b) Folosiți expresia dată la (a) pentru a afla volumul lui (i) un paralelipiped dreptunghiular cu laturile a, b, c (ii) un con circular drept cu înălțimea h și raza bazei R [Sugestie: Calculul este foarte simplu cu conul orientat așa cum se arată în figură ] (iii) o sferă cu raza R Probleme II- (a) Să considerăm o funcție vectorială cu proprietatea V • F = peste tot pe două suprafețe închise S și S și în volumul V cuprins de acestea (vezi figura) Să se arate că fluxul lui F prin Sj este egal cu fluxul lui F prin S În calcularea fluxurilor, alegeți direcția normalelor așa cum este indicată de săgețile din figură (b) Având în vedere câmpul electrostatic al unei sarcini punctiforme q situată la r = , те r r' unde r = x + y + z , arătați prin calcul direct că V • E = , pentru toate r + (c) Demonstrați legea lui Gauss pentru câmpul unei singure sarcini punctiforme date în (b) [Sugestie: este ușor de calculat fluxul lui E pe o sferă centrată la r = ] (d) Cum ați extinde această dovadă pentru a acoperi cazul unei distribuții arbitrare a taxelor? Integrale de suprafață si Divergenţă - (a) Arătați prin calcul direct că teorema divergenței nu este valabilă pentru er F(r,e, Ф) = ^, cu suprafața unei sfere de rază R centrată la origine și V volumul închis De ce eșuează teorema? (b) Verificați prin calcul direct că teorema divergenței este valabilă pentru funcția F din partea (a) când S este suprafața Sx a unei sfere cu raza R plus suprafața S a unei sfere cu raza R , ambele centrate la origine , iar V este volumul cuprins de și S (c) În general, ce restricție trebuie pusă pe o suprafață S, astfel încât teorema de divergență să fie valabilă pentru funcția părții (a)? Capitolul III Integrale de linie și curl A greși de pe calea cea dreaptă este comună omenirii Sofocle Integrale de lucru și linie Am remarcat mai sus că forma diferențială a legii lui Gauss, Ecuațiile ( - ) și ( - ), deși îndeplinește scopul nostru de a lega o proprietate a câmpului electric (divergența acestuia) într-un punct la un punct cunoscut cantitatea (densitatea de sarcină) în același punct, cu toate acestea nu oferă o modalitate convenabilă de a găsi E Motivul este că V • E = p/e este (sau pare să fie) o singură ecuație diferențială cu trei necunoscute (Ел, Ey, Ez) Dar există o altă caracteristică a câmpurilor electrostatice care nu a jucat încă un rol explicit în discuția noastră și care va genera o relație între componentele lui E Astfel, ne va oferi ultimul pas crucial în obținerea unei modalități utile pentru a calcula câmpuri În procesul ex- Integrale de linie și Bucleul Rezolvând această întrebare, vom întâlni unele dintre cele mai importante subiecte în calculul vectorial Proprietatea câmpurilor electrostatice pe care vom începe acum să o discutăm este strâns legată de problema muncii și energiei Vă amintiți fără îndoială definiția elementară a muncii ca forță înmulțită cu distanța Astfel, intr-o dimensiune, daca o forta F(x) actioneaza de la x = a la x = b, munca efectuata este, prin definitie, F(x) Pentru a putea face față unor situații mai generale, trebuie să introducem acum conceptul de integrală de linie Figura III- Să presupunem că avem o curbă C în trei dimensiuni (Figura III- ) și să presupunem că curba este direcționată Prin aceasta ne referim că punem o săgeată pe curbă și spunem "Aceasta este direcția pozitivă" Fie л lungimea arcului măsurată de-a lungul curbei dintr-un punct arbitrar de pe ea cu r = într-un punct P} și s = s la P Să presupunem în continuare că avem o funcție f(x, y, z) definită peste tot pe C Acum să subdivizăm porțiunea lui C între Pt și P arbitrar în N secțiuni Figura III- prezintă un exemplu de astfel de subdiviziune pentru N = Apoi, uniți punctele de subdiviziune succesive prin coarde, dintre care unul tipic, să zicem Zth, are lungimea Asz Acum evaluați/(x, y, z) la (xz, yz, zz), care este orice punct din subdiviziunea Z a curbei și formați produsul /(xz, yz, zz) A yz Făcând acest lucru pentru fiecare dintre cele N segmente ale lui C, formăm suma E/febZ/) Asz Z= Lucru și linie Prin definiție, integrala dreaptă off(x, y, z) de-a lungul curbei C este Integrale limitează această sumă pe măsură ce numărul de subdiviziuni N se apropie în finitate și lungimea fiecărui acord se apropie de zero: f N f(x,y,z)ds = lim zz) CN^"> = fiecare Ss/-" Pentru a evalua integrala de linie, trebuie să cunoaștem calea C De obicei, cea mai convenabilă modalitate de a specifica această cale este parametrică în termeni de parametrul lungimii arcului s Astfel, scriem x = x(s), У - y(s) și z = z(s) Într-o astfel de situație, integrala dreaptă poate fi redusă la o integrală definită obișnuită: [ f(x, jc y, z) ds f[x(s), y(s), z(j)] ds Un exemplu de integrală de linie va fi util aici Pentru simplitate, să lucrăm în două dimensiuni și să evaluăm (x + y) ds, unde C este linia dreaptă de la origine la punctul ale cărui coordonate sunt ( , ) (Figura ІП- ) Dacă (x, y) sunt coordonatele lui Figura III- orice punct P pe C și dacă s este lungimea arcului măsurată de la origine, atunci x = s!\^ și у = sA/ Непсе, x + у = s V = y/ls Prin urmare, + y) ds s ds = V Integrale de linie și Bucleul Figura III- Să integrăm aceeași funcție (x + y) de la ( , ) la ( , ) de-a lungul unei alte căi așa cum se arată în Figura Ш- Nu împărțim integrarea în două părți, una de-a lungul Ct și a doua de-a lungul C Pe Cj avem x = s și у = Astfel, pe C ; x + у = s, și așa De-a lungul C , x = și у = s [rețineți că lungimea arcului pe acest segment al traseului este măsurată din punctul ( , )] Rezultă atunci că f (x + y) ds = [ ( + v) ds = j C J o Adăugând rezultatele pentru cele două segmente, găsim f (x + y) ds = I (x + y) ds + I (x + y) Js = + | = Jc Jc, Jc Lecția de învățat este următoarea: valoarea unei integrale de linie poate depinde (într-adevăr, de obicei este) de calea integrării Integrale de linie care implică funcții vectoriale Deși discuția anterioară ne spune ce este o integrală de linie, tipul de integrală de linie cu care trebuie să ne ocupăm aici are o caracteristică nemenționată încă Vă veți aminti că am introdus discuția noastră despre integralele de linie cu conceptul de muncă Munca, în cea mai ele- Integrale de linie Implicând Vector Funcții simțul mental, este forța ori deplasarea Că aceasta necesită o elaborare devine clar atunci când recunoaștem că atât forța, cât și deplasarea sunt vectori Astfel, luați în considerare o cale C în trei dimensiuni (Figura III- ) Să presupunem că sub acțiunea unei forțe un obiect se deplasează pe această cale de la i, la s În orice punct P al curbei lăsați forța Figura Ш- care acţionează să fie desemnat f(x, y, z) Componenta lui f care funcționează este, prin definiție, numai cea care acționează de-a lungul curbei, adică componenta tangențială Fie t un vector unitar care este tangent la curba la P ' Atunci munca efectuată de forța în deplasarea obiectului de la acesta la s de-a lungul curbei C este W = I f (x, y,z)' t ds, J c unde se înțelege, desigur, că integrarea începe la i = S] și se termină la s = s Noua caracteristică a acestei integrale este că integrând este produsul scalar al două funcții vectoriale Pentru a putea gestiona o astfel de integrală de linie, trebuie să știm cum să găsim t și tocmai la această problemă ne întoarcem acum Considerăm o curbă arbitrară C (Figura III- ) parametrizată după lungimea arcului său La un punct s al curbei avem x = x(s), у = y(s) și z = z(i) În alt punct î + Deoarece avem x + Ax = x(i + Ді), у + Ду = y(s + Ai) și z + Az = z(i + Ai) Prin urmare, t este o funcție a lui x, y și z și ar trebui să fie scris cu adevărat t (x, y, z) Scriem simplu t pentru a evita complicarea notației Integrale de linie și Bucleul Figura III- coarda care unește cele două puncte ale curbei îndreptate de la primul la al doilea este vectorul Ar = iAx + j Ду + kAz, unde Ax = x(s + Ai) - x(i), Ay = y(s + Ai) - y(i), Az = z(i + Ai) - z(i) Dacă acum împărțim acest vector la As, obținem Ar = Topor, Az Ai Ai J Ai Ai' Luând limita acestui lucru pe măsură ce Ai se apropie de zero, rezultă dx, dy | dz ds ds ds' și afirmăm că acesta este t Pentru început, este clar că pe măsură ce Ai -> , vectorul Ar devine tangent la curba de la i Mai departe, în limita Ai -" , vedem că |Ar| -> Ai Prin urmare, în limită mărimea acestei mărimi este Rezultă atunci că putem face identificarea î, dx, dy dz t(i) = i- + i - + k - v ds J ds ds Dacă revenim acum la expresia pentru lucrul W și folosim această formulă pentru t, găsim W = dx dy dz , - + j - + к - ds ds ds ds , de asemenea), obținem lim ~ [fi Д ->o Д J despre (x,y,z) Л dFy F-tds = ~ dx FX mori (III- ) unde $ este notația noastră semicomică, adică circulația în jurul micului dreptunghi S-ar putea să vă întrebați despre generalitatea și unicitatea acestui rezultat, deoarece este obținut folosind o cale de integrare care este specială în două moduri: în primul rând, este un dreptunghi și, în al doilea rând, este paralel cu planul xy Dacă calea nu ar fi un dreptunghi, ci o curbă plană de formă arbitrară, nu ar afecta rezultatul nostru (vezi Problemele III- și III- ) Dar rezultatul nostru depinde cu siguranță de orientarea specială a căii de integrare Alegerea orientării făcută mai sus sugerează în mod clar alte două, iar acestea sunt prezentate în Figura III- (a) și (b) împreună cu rezultatul calculului lim ~ [f| F • t ds ді-"о Д T despre (x,y,z) pentru fiecare Bucul Figura III-l (b) Fiecare dintre aceste trei căi este denumită în onoarea vectorului normal zonei închise Convenția pe care o folosim este următoarea: Trasați curba C astfel încât zona închisă să fie întotdeauna la stânga [Figura III- (a)] Apoi alegeți normalul astfel încât să indice "sus" în di- Figura ПІ-ІЗ(а) Integrale de linie și Bucleul Figura III- (b) recție prezentată în pictare Această convenție este uneori numită regula mâinii drepte, deoarece dacă mâna dreaptă este orientată astfel încât degetele să se îndoaie în direcția în care este trasată curba, degetul mare, întins, arată în direcția normalului [Figura III - (b)J Folosind regula mâinii drepte, avem următoarele: Calcularea ітД ^, $ F • t ds/LS dF, dFy pentru o cale a cărei normală este i, obținem , dFx âFz pentru o cale a cărei normală este i, obținem -r- , r dz dx dFy dFx pentru o cale a cărei normală este k, obținem -r- dx dy (III- a) Rezultă că aceste trei mărimi sunt componentele carteziene ale unui vector Acestui vector îi dăm numele "curl of F", pe care îl scriem curl F Astfel, avem bucla F = i dy ЭГЛ (dFx dFz\ dz ) + \ dz dx J , dFy dFx\ + k - \ dx dy I (III- b) Această expresie este adesea (într-adevăr, de obicei) dată ca definiție a buclei, dar preferăm să o considerăm doar ca forma buclei în coordonate carteziene Vom defini bucla ca limită Bucul Figura III- de circulație către zonă deoarece zona tinde spre zero Pentru a fi precis, fie c"F • t ds circulația lui F pe o cale a cărei normală este n așa cum se arată în Figura III- Apoi, prin definiție n • curl F = lim г!Ѳ! Іпф Prin urmare VXF = sec ф r sin ф g г Ѳ fără apă Sensul buclei Discuția anterioară vă poate lăsa cu sentimentul că a ști cum să definiți și să calculați ondularea unei funcții vectoriale este departe de a ști ce este Faptul că bucla are ceva de-a face cu o linie integrală în jurul unei căi închise (într-adevăr, cuvântul "ondulă" în sine) vă poate sugera că are cumva legătură cu lucrurile care se rotesc, se învârtesc sau se învârtesc Prin intermediul câtorva exemple luate din mișcarea fluidelor, vom încerca să facem puțin mai clare aceste de impresii vagi Înțelesul lui Bucul Figura III- Să presupunem că apa curge pe căi circulare, ceva asemănător cu apa care se scurge dintr-o cadă Un volum mic de apă într-un punct (x, y) la momentul t are coordonatele x = r cos ut, у = r sin wt, unde w este viteza unghiulară constantă a apei (Figura III- ) Astfel , viteza sa la (x, y) este v = i(dx!dt) + \{dyldt) = ro>[-i sin cot + j cos cot] = VXF = Răspuns: Da Exemplul Ar putea fi F = A'f iy - jx) un câmp electrostatic? În acest caz XF = = - k=>VXF = - kK К \ dx dy I Răspuns: Nu Teorema lui Stokes Din aceste exemple putem vedea cât de ușor este de aplicat acest criteriu Teorema lui Stokes Pentru restul acestui capitol, ne îndepărtam de la prezentarea noastră pentru a discuta o altă teoremă faimoasă, una care amintește puternic de teorema divergenței și totuși, după cum vom vedea, destul de diferită de ea Această teoremă, numită după matematicianul Stokes, relaționează o integrală de linie în jurul unui traseu închis la o integrală de suprafață peste ceea ce se numește o suprafață de limitare a căii, așa că primul punct de pe ordinea de zi este definirea acestui termen Să presupunem că avem o curbă închisă C, așa cum se arată în Figura III- (a), și să ne imaginăm că este alcătuită din Figura III- (a) sârmă Acum să presupunem că atașăm o membrană elastică la sârmă, așa cum este indicat în Figura III- (b) Această membrană este o acoperire Figura III- (b) suprafața curbei C Orice altă suprafață care poate fi formată prin întinderea membranei este, de asemenea, o suprafață de acoperire; un exemplu este prezentat în Figura III- (c) Figura III- prezintă patru suprafețe de acoperire diferite ale unui traseu circular plan: (a) regiunea Figura III- (c) Integrale de linie și Bucleul Figura Ш- plan închis de cerc, (b) o emisferă cu cercul drept margine, (c) suprafața curbată a unei calote (un con circular drept) și (d) suprafețele superioare și laterale ale unei cutii de pește de ton Având în vedere aceste remarci preliminare, nu veți fi surprinși să vedem că începem această discuție despre teorema lui Stokes luând în considerare o curbă închisă C și o suprafață de acoperire [Figura III- (a)] Așa cum am făcut anterior, aproximăm această limitare Figura III- (a) suprafață de un poliedru de N fețe, fiecare dintre ele tangentă la S la un punct [Figura III- (b)] Rețineți că acest lucru va crea automat un poligon [marcat cu P în Figura III- (b)] care este de aproximativ Figura III- (b) Teorema lui Stokes Imație la curba C Fie F(x, y, z) o funcție vectorială bine comportată definită în întreaga regiune a spațiului ocupată de curba C și suprafața ei de acoperire S Să formăm circulația lui F în jurul lui С/, limita fetei Z a poliedrului: F • t ds Dacă facem acest lucru pentru fiecare dintre fețele poliedrului și apoi adunăm toate circulațiile, afirmăm că această sumă va fi egală cu circulația lui F în jurul poligonului P: X Ф F • t ds = F • t ds ASh Ci (III- ) unde ASZ este aria feței Z Mărimea dintre paranteze pătrate este, aproximativ, egală cu п/ ( VXF), unde nz este normala pozitivă unitatea de pe fața Z și (VXF)z este curba funcției vectoriale F evaluată în punctul de pe fața Z la care este tangentă la S Spunem "aproximativ" deoarece este de fapt limita ca AS, tinde spre zero a mărimii între paranteze din ecuația (Ш- ), care trebuie identificată ca nz • (VXF)z Ignorând această lipsă de rigoare, scriem N lim l= ] fiecare Д /-> LA FEL DE, F • t ds LA FEL DE, N = lim n, • (VXF), AS; N-+ /-] fiecare Д ,-> (Ш- ) Deoarece curba C este forma limită a poligonului P, avem de asemenea lim V-*co F • t ds = F • t ds (III- ) fiecare Д /-> Teorema lui Stokes Combinând ecuațiile (III- ), (ПІ- ) și (Ш- ), ajungem, în sfârșit, la teorema lui Stokes: F • î ds = (ПІ- ) unde este orice suprafață care limitează curba C Astfel, în cuvinte, teorema lui Stokes spune că integrala de linie a componentei tangențiale a unei funcții vectoriale pe o cale închisă este egală cu integrala de suprafață a componentei normale a curbei acelei funcții integrate peste orice suprafață de acoperire a căii Teorema lui Stokes este valabilă pentru orice funcție vectorială F care este continuă și diferențiabilă și are derivate continue pe C și S Să lucrăm un exemplu Luați F(x, y, z) = iz + jx - kx, cu C cercul cu raza centrat la origine și situat în planul xy și S partea din planul xy închisă de cerc [vezi Figura III- (a)] Acum F • t ds = zdx + x dy - x dz Figura III- (a) Astfel, x dy Până acum am avut întotdeauna curbe parametrizate cu lungimea arcului y În această situație, totuși, traseul C este cel mai ușor parametrizat în termenii unghiului Ѳ prezentat în Figura III- (b) Astfel, scriem C ir cos Ѳ dQ = tt, o (Ш- ) unde folosim x = cos Ѳ și у = sin Ѳ Integrale de linie și Bucleul Figura ІП- (Ь) În continuare calculăm i j к VXF = d/dx d/dy d/dz = j + к Z x Suprafața de acoperire aici este o porțiune a planului xy, astfel încât normala unității în direcția pozitivă este n = k Prin urmare, n • VXF = к • ( j + к) = și • VXF dS = (III- ) unde această ultimă egalitate rezultă din faptul că integrala de suprafață în acest caz este doar aria cercului unitar Deoarece acest rezultat [Ecuația (III- )] este identic cu cel obținut mai sus [Ecuația (III- )], am ilustrat teorema lui Stokes [Ecuația (III- )] Să refacem acest calcul, de data aceasta alegând emisfera prezentată în Figura III- (c) ca suprafață de acoperire S Folosind ecuația ( - ) cu F înlocuit cu VXF, obținem + dx dy = j J ~dxdy + jj dxdy O aplicație a teoremei lui Stokes z Figura Ш- (с) unde R este cercul unitar în planul xy prezentat în Figura III- (a) A doua integrală din cea mai dreaptă egalitate de mai sus este doar aria acelui cerc, deci valoarea sa este ir Prima integrală poate fi gestionată prin introducerea coordonatelor polare Găsim dx dy Este ușor de arătat că integrala peste Ѳ este zero Prin urmare, JLfi-v XF dS = tt, în concordanță cu rezultatul nostru anterior O aplicație a teoremei lui Stokes O aplicație importantă a teoremei lui Stokes este furnizată de legea circuitală a lui Am-pere Luați în considerare orice buclă închisă C care include un curent I ca în Figura III- Rețineți că direcția lui C și cea a lui I corespund aceleiași reguli din dreapta care raportează direcțiile lui C și normala pozitivă la o suprafață care limitează C Ampere Figura III- Integrale de linie și Bucleul Legea circuitului spune că integrala de linie a câmpului magnetic В este legată de curent astfel: В • t ds unde constanta po, numită permeabilitatea spațiului liber, are valoarea , X newtoni pe amper Această lege, la fel ca legea lui Gauss și legea circulației, spune ceva despre integrala unui câmp (în acest caz, câmpul magnetic) și, la fel ca în cele două cazuri precedente, este convenabil să o reexprimați astfel încât să fie ne teii ceva despre domeniu la un punct În acest scop, introducem mai întâi densitatea de curent J (vezi pagina ) Astfel, dacă curentul curge printr-o zonă AS cu normal n (Figura III- ), densitatea de curent J este astfel încât M = J • n AS, Figura III- unde AZ este curentul total Adică, densitatea de curent este o funcție vectorială a cărei mărime este curentul pe unitate de suprafață și a cărei direcție este cea a fluxului de curent Dacă J(x, y, z) este densitatea curentului, atunci curentul total care curge printr-o suprafață S este J • n dS Astfel, legea lui Ampere poate fi scrisă В • t ds Ho S poate fi orice suprafață care acoperă curba C Dacă, așa cum este de obicei, curentul trece printr-un fir a cărui secțiune transversală Teorema lui Stokes și regiunile pur și simplu conectate Figura III- nu include întreaga suprafață de acoperire, nu contează; putem integra peste mai mult decât secțiunea transversală a sârmei dacă ne amintim că J pentru acea parte a suprafeței S tăiată de sârmă și J = pentru restul (Figura III- ) Prin urmare, J/j âdS secțiunea transversală a firului ff J"n ds întreaga suprafață de acoperire S Acum folosind teorema lui Stokes [Ecuația (III- )], avem В • t dS = n • J dS Deoarece C și S sunt arbitrare, concluzionăm că VX В = p-qJ Aceasta este forma diferențială a legii lui Ampere Este, de asemenea, un caz special al uneia dintre ecuațiile lui Maxwell, valabilă atunci când câmpurile nu variază în timp Teorema lui Stokes și regiunile pur și simplu conectate Pentru multe scopuri, inclusiv pentru unele aplicații importante, trebuie să putem afirma că teorema lui Stokes este valabilă în întreaga regiune D din spațiul tridimensional Prin aceasta, vrem ca teorema să fie valabilă pentru orice curbă închisă C situată în întregime în D și pentru orice suprafață de acoperire a lui C aflată, de asemenea, în întregime în D Aceasta, desigur, înseamnă că funcția F trebuie să fie continuă și diferențiabilă și să aibă prime derivate continue în D Dar în plus trebuie să impunem o restricție asupra regiunii D însăși A întelege Integrale de linie și Bucleul cum se întâmplă acest lucru, să presupunem mai întâi că D este interiorul unei sfere Dacă F este neted peste tot în D, atunci teorema lui Stokes este valabilă pentru orice curbă închisă C situată în întregime în D și pentru orice suprafață de acoperire a lui C aflată, de asemenea, în întregime în D Cu alte cuvinte, teorema lui Stokes este valabilă peste tot în D Un pic de gândire ar trebui să vă convingă că aceeași linie de raționament se aplică regiunii dintre două sfere concentrice, cu condiția ca F să fie netedă în acea regiune Dar pentru anumite tipuri de regiuni pot apărea probleme Ca exemplu, să presupunem că D este interiorul unui tor (aproximativ ca un bagel sau un tub interior umflat; vezi Figura III- ) Problema în acest caz este că este Figura ІП- este posibil să se construiască o curbă închisă în D ca cea prezentată în figură, cu proprietatea că nici una dintre suprafețele sale de acoperire nu se află în întregime în D Deși insistăm ca F să fie netedă în D, nu i se impune condiții în altă parte, deci că în afara regiunii s-ar putea să nu îndeplinească cerinţele de netezime care asigură va-liditatea teoremei lui Stokes Relația dintre integrala de linie peste C și integrala de suprafață peste afirmată de teoremă poate, și în multe cazuri o face, să se rupă dacă F nu este netedă pe Matematicienii se referă la regiuni precum interiorul unei sfere sau spațiul dintre două sfere concentrice ca fiind pur și simplu conectate, în timp ce interiorul unui tor nu este pur și simplu conectat Prin definiție, o regiune D este pur și simplu conectată dacă orice curbă închisă care se află în întregime în D se poate micșora până la un punct fără a lăsa D Folosind această definiție, ar trebui să puteți verifica dacă interiorul unei sfere și regiunea dintre două sferele concentrice sunt ambele pur și simplu conectate, dar interiorul unui tor nu este Având la dispoziție conceptul de conexiune simplă, putem specifica cu ușurință condițiile în care teorema lui Stokes este valabilă într-o regiune: funcția vectorială F trebuie În continuare, când spunem că o funcție este "netedă", vom spune că este continuă, diferențiabilă și are derivate prime continue Independența căii și bucla să fie netedă peste tot într-o regiune simplu conectată D Atunci teorema lui Stokes [Ecuația (III- )] este valabilă pentru orice curbă închisă C și orice suprafață de acoperire a lui C, ambele fiind în întregime în D De cele mai multe ori vom presupune că funcțiile cu care lucrăm sunt netede și că regiunile de interes sunt pur și simplu conectate Există însă situații, precum cea discutată în secțiunea următoare, în care conexiunea simplă joacă un rol esențial și le vom evidenția pe măsură ce ajungem la ele Independența căii și bucla În discuția noastră despre forma diferențială a legii circulației, am arătat că, deoarece integrala de linie a unui câmp electrostatic E este zero pe orice cale închisă, curba lui E este zero Același lucru este valabil și pentru orice funcție vectorială F; adică dacă F • t ds = pentru toate căile închise C, atunci VXF = Dovada acestui fapt este exact aceeași cu cea dată la paginile - cu E înlocuit peste tot cu F Este adevărat și inversul acestei afirmații? Adică, dacă VXF = , aceasta înseamnă că circulația lui F este zero pe toate căile închise? La prima vedere ar putea părea că răspunsul la această întrebare este da Tot ce trebuie să facem este să folosim teorema lui Stokes și să observăm că, deoarece prin ipoteza VXF = , F • t ds = = s Cu toate acestea, există un defect în această linie de raționament Reamintim că validitatea teoremei lui Stokes necesită ca F să fie netedă într-o regiune pur și simplu conectată Dacă regiunea nu este pur și simplu conectată, teorema lui Stokes s-ar putea să nu fie valabilă, cel puțin pentru unele căi închise care se află în regiune, iar faptul că VXF = nu garantează că circulația lui F este zero pe toate căile închise Cel mai aproape ne putem ajunge la o conversă este să spunem că dacă VXF = peste tot în а Integrale de linie și Bucleul regiune pur și simplu conectată, atunci circulația lui F este zero pentru toate căile închise din acea regiune Cele două afirmații "circulația este egală cu zero" și "ondulare este egală cu zero" sunt echivalente doar într-o regiune pur și simplu conectată Există o modalitate ușor diferită, dar adesea utilă, de a afirma această legătură între circulație și bucle; și anume, dacă fc F • t ds este independent de cale, atunci VXF = , iar dacă VXF = într-o regiune pur și simplu conectată, atunci fc F • t ds este independent de cale Nu ar trebui să aveți dificultăți în a stabili acest lucru pentru dvs PROBLEME Ш-l Utilizați un argument ca cel dat în text pentru forța Coulomb (paginile - ) pentru a arăta că fc F • t ds este independent de calea oricărei forțe centrale F Ш- În text am obținut rezultatul prin integrarea pe un mic drum dreptunghiular Ca un exemplu al faptului că acest rezultat este independent de cale, rederive-l, folosind calea triunghiulară prezentată în figură III- Calculați curba fiecăreia dintre următoarele funcții folosind ecuația (III- b): (a) iz + jx - ky (b) ixz - kx (c) ie~y + je"z + ke Jt (d) iyz + jxz + kxy (e) -iyz + jxz (f) ix + jy + k(x + y ) (g) ixy + jy + kyz (h) (ix + jy + kz)/(x + y + z ) / , (x, y, z) ( , , ) Probleme III^l (a) Calculați f F • t ds pentru funcția din problema III- (a) pe un drum pătrat de latură centrat pe (x , y , ), situat în planul xy și orientat astfel că fiecare latură este paralelă cu axa x sau y (b) Împărțiți rezultatul părții (a) la aria pătratului și luați limita coeficientului ca s -" Comparați rezultatul cu componenta z a buclei găsite în problema III- (a) (c) Repetați părțile (a) și (b) pentru funcțiile din Problema III- (b), (c) și (d) (Puteți găsi interesant să încercați căi de diferite orientări și/sau forme ) - (a) Calculați $ F • t ds unde F = k(y + y ) peste perimetrul triunghiului prezentat în figură (integrați în direcția indicată de săgeți) (b) Împărțiți rezultatul părții (a) la aria triunghiului și luați limita ca a -" (c) Arătați că rezultatul părții (b) este n • VXF evaluat la ( , , ) unde fi este vectorul unitar normal triunghiului și îndreptat departe de origine Ш- Arată că VXA Г = A, unde r = ix + jy + kz și A este un vector constant III- Să se arate că V • (V x F) = (Să presupunem că derivatele parțiale secundare mixte sunt independente de ordinea diferențierii De exemplu, d F Jdx dz = d F /dz Эх ) Ш- În text (paginile - ) am obținut componenta z a lui VXF în coordonate cilindrice Procedând în același mod, obțineți componentele Ѳ și r date la pagina - Urmând procedura sugerată în text (paginile - ), obțineți expresia pentru VXF în coordonate sferice prezentate la pagina Cifrele date la pagina vă vor fi de ajutor Integrale de linie și Bucleul ПІ- (a) Rescrieți funcția din problema III- (e) în coordonate cilindrice și calculați curba acesteia folosind expresia dată la pagina Convertiți rezultatul înapoi în coordonate carteziene și comparați cu răspunsul obținut în problema III- (e) (vezi problema П- ) (b) Repetați calculul de mai sus pentru funcția Problemei ID- (f) III-ll (a) Rescrieți funcția din problema III- (g) în coordonate sferice și calculați curba acesteia folosind expresia dată la pagina Convertiți rezultatul înapoi în coordonate carteziene și comparați cu răspunsul obținut în problema HI- (g) (vezi problema П- ) (b) Repetați calculul de mai sus pentru funcția Problemei III- (h) z Ш- Orice forță centrală poate fi scrisă sub formă F(r) = ej(r), unde er este un vector unitar pe direcția radială și/este o funcție scalară Arătați prin calculul direct al buclei că această funcție este irrotațională (adică VXF = ) Probleme III- Care dintre funcțiile din problema ПІ- ar putea fi câmpuri electrostatice? III- Folosiți teorema lui Stokes pentru a arăta asta t ds = , unde C este o curbă închisă și t este un vector unitar tangent la curba C III- Verificați teorema lui Stokes F • t ds = JJn-V X FdS în fiecare dintre cazurile următoare: (a) F = iz - jy C, pătratul laturii situat în planul xz și direcționat așa cum se arată , cele cinci pătrate Sb S , S , S și S așa cum se arată în figură (b) F = iy + jz + kx C, cele trei sferturi de cerc Cb C și C direcționate așa cum se arată în figură S, octantul sferei x + y + z = înconjurat de cele trei arce X Integrale de linie și Bucleul (c) F = iy - jx + kz C, cercul cu raza R situat în planul xy, centrat pe ( , , ) și direcționat așa cum se arată în figură S, suprafețele curbe și superioare ale cilindrului cu raza R și înălțimea h - (a) Să considerăm o funcție vectorială cu proprietatea VXF = peste tot pe două curbe închise C} și C și pe orice suprafață de acoperire S a regiunii cuprinse de acestea (vezi figura) Să se arate că circulația lui F în jurul lui C{ este egală cu circulația lui F în jurul lui C În calcularea circulațiilor, direcționați curbele așa cum este indicat de săgețile din figură (b) Câmpul magnetic datorat unui fir drept infinit lung care transportă un curent uniform I este В = (р оІ/ тгг)ев Arătați că VX В = peste tot, cu excepția r = (c) Demonstrați legea circuitului lui Ampere pentru câmpul firului dat în partea (b) [Sugestie: Folosiți rezultatul (b) pentru a găsi circulația lui В Probleme în jurul unui cerc cu firul care trece prin centru și normal cu planul său Apoi utilizați rezultatul din partea (a) pentru a reimite această circulație la circulația în jurul unei curbe arbitrare care înconjoară curentul ] ІП- (a) Se consideră funcția dată în coordonate cilindrice de ёѳ F(r, Ѳ, z) = -y Arătați că teorema lui Stokes nu este valabilă pentru această funcție dacă C este cercul cu raza R în planul xy centrat la origine, iar este porțiunea planului xy închisă de C De ce eșuează teorema în acest caz? (b) Luați în considerare regiunea D care constă din tot spațiul tridimensional cu axa z îndepărtată Este funcția F definită în (a) netedă în D? Teorema lui Stokes este valabilă în D? Este D o regiune pur și simplu conectată? III- Forța electromotoare 'в într-un circuit C este egală cu circulația câmpului electric E în jurul circuitului: g = E • t ds ' c Faraday a descoperit că într-un circuit staționar o forță electromotoare este indusă de un flux magnetic în schimbare Acesta este, dt Unde Ф = JJ В • n dS, t este timpul (nu îl confundați cu vectorul tangent t) și S este orice suprafață de limitare a lui C Folosiți această informație și teorema lui Stokes pentru a deriva ecuația VXE = dB dt ' care este una dintre ecuațiile lui Maxwell Ш- Să se determine valoarea integralei drepte fc F • t ds, unde F = (e~y - ze x)i + (e~z - xe~-')j + (e~x - ye z)k iar C este calea % = hT ln( +p)' Ttp У = păcat -, F • t ds = X Ф F ' t ds p = •> c = C(x , y ) ДА + (tiv) ix,-x ) SUA, \C'X/xoJ"(=l + (¥) (Л-Уо)Д , + - \оу/х",л/=і unde C, este perimetrul zonei Z, (x , y ) este un punct din regiunea închisă de P și ДА este aria cuprinsă de P (c) Arată că lim N-*°o fiecare AS,-> > F • t ds = Ф F • t ds p J c C(x , y ) + (x - x )l unde Д este aria regiunii R închisă de C și (x, y) sunt coordonatele centroidului regiunii R; acesta este, x = WJRxdxdy și y = XsllRydxdy- (d) În cele din urmă, calculați (VXF), = lim - x dy dz dy\dzj dz\dy J Э ф Э ф = - = dy dz dz dy Această ultimă egalitate urmează dacă ф și derivatele sale prima și a doua sunt continue, pentru că atunci д ф/ду dz = д ф/дг dy Evident, se poate demonstra că celelalte două componente ale VXF dispar exact în același mod Prin urmare, Эф = ~ (q=x,y,z) => VXF = Reversul a ceea ce tocmai am arătat ar afirma că dacă VXF = , atunci există o funcție scalară ф astfel încât F = Ѵф, o afirmație care este adevărată, cu condiția ca regiunea de interes să fie pur și simplu conectată Pentru a înțelege acest lucru, putem consulta Figura IV- , care arată cum independența de cale a integralei drepte a lui F • t, VXF = și F = Ѵф sunt legate Săgețile solide din diagramă reprezintă implicații care sunt valabile în general, cu condiția ca F să fie neted Săgețile întrerupte reprezintă implicații care necesită nu numai ca F să fie netedă, ci și ca regiunea de interes să fie pur și simplu conectată Am arătat deja că ( ) implică atât ( ) cât și ( ) și că ( ) implică ( ) într-o regiune pur și simplu conectată Combinând aceste două afirmații, vedem că ( ) implică ( ) într-o regiune pur și simplu conectată În practică, la fel cum funcțiile cu care ne ocupăm de obicei au derivate prime continue (și, prin urmare, sunt netede), regiunile cu care lucrăm sunt pur și simplu conectate În astfel de circumstanțe, ne putem relaxa puțin și luăm în considerare cele trei afirmații rezumate în fig Gradientul independent de cale VxF = Figura IV- ure IV- ca echivalent: fiecare implică și este subînțeles de fiecare dintre celelalte Cu toate acestea, ar trebui să fiți conștienți de conexiunea simplă și de implicațiile acesteia pentru relațiile dintre cele trei afirmații Pentru a da un exemplu simplu al ideilor pe care le-am discutat, luăm în considerare funcția vectorială F(x, y, z) = iy + jx Această funcție este netedă peste tot și am observat deja că curba ei este zero (pagina ) Conform a ceea ce tocmai am spus, aceasta înseamnă că trebuie să existe o funcție scalară ф(х, у, z) astfel încât F să fie gradientul ei Astfel, ф trebuie să satisfacă p Эф Эф F = = - г dz' În mod clar, ф(х, у, z) = xy + C, unde C este o constantă arbitrară, satisface aceste relații Acest lucru ar trebui să fie în contrast cu cazul funcției F = iy - jx, a cărei ondulare nu dispare (pagina ) Dacă această funcție ar fi gradientul unei funcții scalare ф, ar trebui să avem Эф Эх ' P дф Fy=~x = ~dr Эф F = = , z dz Găsirea electrostaticului Camp dar, așa cum ar trebui să vă puteți convinge, nu există nicio funcție ф care să satisfacă aceste trei ecuații Expresia pe care am scris-o pentru gradientul unei funcții scalare ф(х, у, z), și anume, Ѵф = дф Эф , Эф -J -к -, dx dy dz este de fapt doar forma acestui operator în coordonate carteziene Pentru a găsi forma gradientului în alte sisteme de coordonate, dacă o procedați direct, este o muncă plictisitoare De exemplu, pentru a găsi gradientul în coordonate cilindrice, ar trebui mai întâi să exprimăm vectorii unitari cartezieni i, j și к în termeni de mărimi analoge er, ёѳ și e, în coordonate cilindrice Atunci, folosind x = r cos Ѳ, у = r sin Ѳ, și regula lanțului pentru diferențiere, ar trebui să exprimăm derivatele față de x, y și z în termenii celor față de г, Ѳ, și z- Nu vom urmări această chestiune aici, deoarece mai târziu (vezi paginile și urm ) ne va fi disponibilă o metodă mai simplă și mai rapidă În prezent, menționăm doar forma gradientului în coordonate cilindrice și sferice Cilindric: Эф g Эф r dr еѳ г ЭѲ л Эф вг dz' (IV- ) Sferic: л Эф Эф л Эф ѵф - е, -г-I- вфу-г + еѳ-;-тзз дг ф г Эф г sin ф ЭѲ (ІѴ^) O definiție fără coordonate a gradientului analogă celor date pentru divergența [Ecuația ( - )] și curl [Ecuația (III- )] este discutată în Problema IV- Găsirea câmpului electrostatic Am început discuția despre calculul vectorial cu o căutare a unei metode convenabile pentru găsirea câmpului electrostatic Investigațiile noastre ne-au condus la forma diferențială a legii lui Gauss, V • E = p/e Gradientul Chiar și această expresie nu este adesea utilă pentru a găsi E, deoarece este o ecuație în trei necunoscute (Ex, Ey și E în coordonate carteziene) Acum, în sfârșit, suntem capabili să ne încheiem discuția și să scriem ecuațiile care sunt adesea cele mai utile dintre toate metodele cunoscute pentru găsirea câmpului Acest ultim pas se bazează pe observația că din moment ce E • t ds = pentru orice cale închisă C, câmpul E poate fi scris ca gradient al unei funcții scalare În mod convențional, această funcție, numită potențial electrostatic, este desemnată Ф(х, у, z), și scriem E = -ѴФ Combinând această ecuație cu forma diferențială a legii lui Gauss [Ecuația ( - )], obținem Ѵ'(-ѴФ) = p/e , sau V • (ѴФ) = - p/e Când scriem partea stângă a acestei ecuații în detaliu, găsim Ѵ-(ѴФ)= i^- + j^- + k|-\ dx dy dz дФ , ЭФ , ЭФ i + J , + к - dx dy dz Э Ф + Э Ф /dx\ și astfel ecuația lui Laplace și condițiile la limită asociate sunt Această secțiune nu este esențială pentru ceea ce urmează și poate fi omisă Folosind a lui Laplace Ecuaţie Figura IV- și O Vo atx = O la x = s Aceasta este o problemă banală și soluția este Câmpul electric se găsește folosind E = -ѴФ, care dă Ex = -și Ey = E, = Astfel, câmpul este un vector constant normal plăcilor Aceasta este o excelentă aproximare a potențialului și a câmpului dintre, dar departe de marginile, două plăci ale căror dimensiuni liniare sunt mari în comparație cu separarea lor Puteți recunoaște acest aranjament ca un condensator cu plăci paralele Al doilea exemplu este un condensator sferic, adică două sfere concentrice cu raze Rt și R , cu cea interioară menținută la o potențială Vo și cea exterioară la zero (Figura ГѴ- ) Ni se cere să găsim potențialul și câmpul peste tot între sfere În această situație, evident, am face bine să muncim Gradientul Figura IV- în coordonatele sferice г, Ѳ și ф, în care ecuația lui Laplace dintre sfere are forma impunătoare Ѵ Ф = іа/ \ Ф(г, Ѳ) = -Eor - cos Ѳ Pentru a găsi câmpul electric, procedăm ca de obicei cu E = - ѴФ Folosind ecuația (IV- ), obținem p - -дФ - F Er- ~^-E R\ cos Ѳ, ц ЭФ r г ЭѲ păcat Ѳ, - R E = = z dz Ar trebui să verificați că pentru r mare, acest câmp se reduce la Eoi după cum este necesar S-ar putea să găsiți acest ultim exemplu neliniştitor, deoarece o anumită cantitate de ghicituri profunde este folosită pentru a găsi potențialul De fapt, există proceduri standard care, în probleme de acest gen, duc mai mult sau mai puțin direct la rezolvare O discuție despre acestea Direcțional Derivatives and the Gradient procedurile, totuși, ar fi foarte lungi și (în expresia de bine) depășesc domeniul de aplicare al acestui text Înainte de a trece mai departe, totuși, merită menționat un alt punct: o soluție a ecuației lui Laplace care satisface condițiile la limită adecvate este unică Adică, există una și o singură astfel de soluție, astfel încât, dacă rezolvăm o problemă prin ghicituri și skulluggery, iar altcineva o rezolvă cu tehnici matematice rafinate și elegante, cele două Soluții, în ciuda pedigree-urilor lor disparate, trebuie să fi la fel În Problema IV- veți fi conduși printr-o dovadă a acestui fapt remarcabil Direcțional Derivatives and the Gradient Am introdus gradientul ca un fel de artificiu matematic util în discutarea integralelor de linii independente de cale Ne întoarcem acum la o examinare mai detaliată a gradientului pentru a descrie semnificația sa geometrică Înainte de a începe discuția noastră, facem câteva comentarii despre seria Taylor, deoarece acestea sunt necesare în cele ce urmează Pentru o funcție scalară a unei variabile care este adecvat continuă și diferențiabilă, avem f(x + Дх) = /(x) + Дх/'(х) + |(Дх) /"(х) + ■ ■ • Aceasta spune că valoarea funcției la un anumit punct x + Дх poate fi scrisă ca suma (de obicei) a infinitului de termeni care implică funcția și derivatele ei în alt punct x Printre altele, această serie Taylor este utilă pentru calcul, deoarece dacă cele două puncte sunt apropiate (adică dacă Дх este mic), atunci putem trunchia seria după un anumit număr de termeni (care sperăm că este mic), întrucât termenii neglijați, fiecare proporțional cu o putere mare a numărului mic Дх, se vor însuma la o valoare care este neglijabilă Seria Taylor poate fi formată și pentru funcții ale mai multor variabile Astfel, pentru o funcție a două variabile avem /(x +Дх, y + Ду) = /(х,у) + Дх|+Ду| + - (IV- ) Aceasta spune că valoarea funcției la un moment dat (x + Дх, у + Ду) poate fi scrisă ca o sumă a (de obicei) infiniti termeni care implică funcția și derivatele ei într-un alt punct Gradientul (X y) Nu vom avea niciodată nevoie de forma explicită a termenilor rămași din această serie [reprezentată prin puncte în ecuația (IV- )] Ar trebui să știm, totuși, că acești termeni implică puteri mai mari ale numerelor "mici" Ax și Ay (de exemplu, Ax , Ay , AxAy, Ax , Ay , Ax Ay și așa mai departe) Având în vedere aceste idei simple, ne întoarcem acum la sarcina noastră principală Să considerăm o funcție z = fix, y) Din punct de vedere geometric, aceasta reprezintă o suprafață așa cum se arată în Figura IV- (a) Fie (x, y) coordonatele unui punct P din planul xy Inaltimea suprafetei deasupra acestui punct este reprezentata de lungimea liniei punctate PQ-, adica PQ = z~ fix, y) Să presupunem că acum facem un pas scurt în planul xy către un nou punct P' cu coordonate (x + Ax, у + Ay) Înălțimea suprafeței deasupra acestui punct este P'Q' = fix + Ax, у + Ay) Fie Av lungimea pasului (Ax = PP') Figura IV- (a) În continuare, întrebăm cât de mult s-a schimbat funcția/s-a schimbat ca urmare a efectuării acestui pas În mod clar, această schimbare este diferența dintre cele două înălțimi PQ și P'Q', și P'Q' ~ Pg = Д/ = /(x + Ax, у + Ay) - fix, y) Aplicând formula seriei Taylor menționată mai sus [Ecuația (IV- )], obținem ^ +м +- dx J dy Direcțional Derivatives and the Gradient Acum reformulam această expresie prin ceea ce la început poate părea o elaborare inutilă a notației Fie As un vector care are mărimea Ai și puncte de la P la P' Clar, As = i Ax + j Ay Dar gradientul de/este dx + J dy (o specializare evidentă a notației gradient la o funcție a două, mai degrabă decât trei, variabile) Urmează deodată că Af= (As) •(?/) + •••■ Complicarea contează puțin mai mult, să fie u un vector unitar în direcția lui As Apoi As = u Ai și A/= (u • V/) Ai + •■ ■, astfel încât u"V/+ - Ai J Acum luăm limita acestei ecuații pentru a obține ds hm - = u • Vf Ал- Ai (IV- ) Nu mai este nevoie de "+ • ••", deoarece punctele reprezintă termeni care merg la zero pe măsură ce Ai merge la zero Această nouă expresie [Ecuația (IV- )] are o interpretare simplă: este rata de schimbare a funcției/(x, y) în direcția lui As (adică a lui u) Redesenând figura IV- (a) și trecând un plan prin P și P' paralel cu axa z [Figura IV- (b)], vedem că decupează suprafața z - f(x, y) într-un curba C Mărimea df/ds definită în ecuația (IV- ) este panta acestei curbe în punctul Q Gradientul Figura ГѴ- (Ь) Mărimea df/ds se numește derivată directă a lui f Deși analiza dată mai devreme care a condus la această derivată a fost pentru funcții a două variabile, toate rezultatele se aplică funcțiilor a trei (sau mai multe) variabile Prin urmare, ^F(x,y,z) = u*VF la fel de este rata de schimbare a funcției Fix, y, z) în direcția specificată de vectorul unitar u Un exemplu de derivat direct poate fi amuzant aici Vom lucra cu o funcție a două variabile, astfel încât să putem desena imagini Deci, să luăm în considerare z=f(x,y) = (x +y ) ' , care este un con circular drept inversat a cărui axă coincide cu axa z [vezi Figura IV- (a)] Solicităm derivata directă a acestei funcții la un punct x = a și у = b și în direcția specificată de u = i cos Ѳ + j sin Ѳ [vezi Figura IV- (b)] Mai întâi avem nevoie de gradientul fix, y) Dar У X , df У дх Z și dy z ' după cum puteți verifica cu ușurință Prin urmare, " ix + iy Direcțional Derivatives and the Gradient Figura IV- (a) Figura IV- (b) și df nr xc°s Ѳ + у sin Ѳ a cos Ѳ + b sin Ѳ = " v/= Să presupunem că Ѳ este ales astfel încât u să fie în direcția radială, așa cum este indicat în Figura IV- (c) Aceasta înseamnă cos Ѳ = A (a + b ) ' ' sin Ѳ = b (a + b ) ' și așa df a a b b ds \/ a + b Va + Z> Va + Z> Va + fr Semnificația acestui rezultat este evidențiată în Figura IV- (d) Gradientul Figura IV- (d) Un al doilea caz interesant este cel în care u este ales perpendicular pe direcția exemplului anterior [vezi Figura IV- (e)] Avem atunci cos Ѳ = - b (a + b )m' sin Ѳ A (a + b ) ' și așa df a / b \ ds Va + b \ Va + b ) Semnificația acestui rezultat este ilustrată în Figura IV- (f) Geometric Semnificația gradientului Figura IV- (e) Figura IV- (f) Semnificația geometrică a gradientului Având la dispoziție conceptul de derivată direcțională, suntem acum în măsură să oferim o interpretare geometrică a gradientului La un moment dat Po cu coordonatele (x , y , z ) avem unde indicele înseamnă că cantitatea trebuie evaluată în punctul (x , y , to) Acum (VF) , gradientul lui F evaluat la Po, poate fi reprezentat printr-o săgeată care emană din acel punct, așa cum se arată în Figura IV- (a) Dacă întrebăm în ce direcție trebuie să ne mișcăm pentru a face (dF/ds)a cât mai mare posibil, este clar că u ar trebui să fie în aceeași direcție cu (VF) Acest lucru se datorează faptului că dacă lăsăm a fi unghiul dintre u și (VF) , atunci (dF!ds)Q = |VF| cos a, și acesta este cât de mare poate fi atunci când a = Astfel, gradientul de dj /o u' (W Gradientul Figura IV- (a) o funcție scalară F(x, y, z) este un vector care se află în direcția în care F suferă cea mai mare rată de creștere și care are o mărime egală cu rata de creștere în acea direcție Pentru a ilustra această interpretare a gradientului, să ne întoarcem la conul inversat z = /(x, y) = (x + y ) / despre care am discutat mai devreme Am învățat asta V/ = к + jy z și ds a cos Ѳ + & sin Ѳ Va + b = ( ) Pentru a găsi direcția în care/(x, y) suferă cea mai mare rată de schimbare, setăm dD -a sin + cos " + b Aceasta dă tan = bl a, de unde cos = a/(a + b )' și sin = bl (a + b ) Deci (dflds)mm = Pe de altă parte, |V/I = = , deoarece z = x + y În plus, tan = b/a corespunde direcției ai + bj, în timp ce în punctul (a, b), W = ai + bj (a + b ) ' ' Geometric Semnificația gradientului care este un vector în aceeași direcție Astfel, ambele proprietăți ale gradientului sunt ilustrate; este în direcția ratei maxime de creștere, iar magnitudinea sa este egală cu rata de creștere în acea direcție Ca un al doilea exemplu de interpretare a gradientului, considerăm planul z = fix, у) = - x - у prezentat în Figura IV- (b) Este ușor de observat că V/= -i - j Folosind u = i cos Ѳ + j sin Ѳ ca înainte, se tind df/ds = u • V/= -cos Ѳ - sin Ѳ = ( ) Prin urmare ~~ = sin - cos = u Figura IV- (b) care dă O = tt/ sau tt/ Testul derivatei a doua arată că ir/ corespunde unui minim și tt/ , unui maxim Din figura IV- (b) se poate observa că cea mai mare rată de creștere este într-adevăr la un unghi de tt/ Mai mult, cea mai mare rată de creștere este = О( тт/ ) = V , întrucât |V/| = |-i - j| = V Din nou, ambele proprietăți ale gradientului sunt ilustrate de acest exemplu Având la dispoziție această interpretare geometrică a gradientului, acum putem vedea motivul semnului negativ în ecuația E = - ѴФ: Deoarece ѴФ este un vector în direcția creșterii Ф, forța asupra unei sarcini pozitive q este F = б/Е = - с/ѴФ, care este în direcția descreșterii Ф Astfel, semnul negativ asigură că o sarcină pozitivă se deplasează "în jos" de la o potențială mai mare la una mai mică Există o altă proprietate a gradientului utilă în înțelegerea semnificației sale geometrice Pentru a concretiza această discuție, să fie T(x, y, z) o funcție scalară care dă temperatura la Gradientul Figura IV- orice punct (x, у, z) Locul tuturor punctelor cu aceeași temperatură To este (în cel mai simplu caz) o suprafață a cărei ecuație este T(x, y, z) = Tb (Figura IV- ) Aceasta se numește suprafață izotermă Acum arătăm că VT este un vector normal la suprafața izotermă Fie C orice curbă situată pe suprafața izotermă și fie P orice punct pe C Fie u vectorul unitar tangent la C la P (nu contează ce direcție de-a lungul C luăm) Derivata direcțională în direcția u este ds) = u • V = deoarece T nu se modifică pe măsură ce ne deplasăm de-a lungul suprafeței izoterme Dacă produsul scalar a doi vectori, niciunul dintre ei zero, dispare, cei doi vectori sunt perpendiculari Astfel, VT este perpendicular pe C la P Prin același argument este perpendicular pe orice curbă de pe suprafața prin P (cum ar fi C în Figura IV- ) Dar acest lucru poate fi adevărat numai dacă VT este normală la suprafața izotermă la P În general, atunci, V/(x, y, z), unde f(x, y, z) este o funcție scalară, este normală la suprafață f(x, y, z) = constant Un exemplu simplu al acestei proprietăți a gradientului este oferit de funcția F(x, y, z) - x + y + z - Suprafața F(x,y,z) = constantă este, desigur, o sferă (presupunând că constanta este pozitivă) După cum ar trebui să verificați singur, VF = (ix + jy + kz) = r Astfel, avem un rezultat familiar: un vector normal la o suprafață sferică este în direcția radială Vă lăsăm să vă gândiți la relația geometrică dintre câmpul electrostatic E și suprafețele sale echipotențiale Ф(х, у, z) = constantă Legătura dintre această proprietate a gradientului și expresia noastră anterioară pentru vectorul unitar normal la o suprafață [Ecuația (II- )] este subiectul problemei IV- Gradientul în Cilindrică și sferice Coordonatele Putem face o legătură simplă între proprietatea gradientului tocmai discutat și faptul că este în direcția celei mai mari rate de creștere Orice deplasare față de fix de suprafață, y, z) = constantă, considerată ca un vector s, poate fi rezolvată într-o componentă de-a lungul suprafeței ($ц) și una normală cu aceasta O'±), așa cum se arată în Figura IV- Acea parte a deplasării de-a lungul suprafeței este "mișcare irosită" dacă scopul nostru în mișcare este de a provoca o schimbare a valorii lui/(x, y, z) Doar componenta normală ne duce departe de suprafață și provoacă o modificare a/ Din aceasta este clar că cea mai mare creștere posibilă pentru o anumită mărime a deplasării ar trebui să aibă loc atunci când ne îndepărtăm de suprafață în direcția normală Dar am stabilit deja că cea mai mare rată de creștere are loc în direcția gradientului Astfel, gradientul este normal la suprafață Figura ГѴ- Gradientul în coordonate cilindrice și sferice Un produs secundar al discuției noastre despre derivata direcțională este metoda "mai ușoară și mai rapidă" pentru calcularea gradientului în coordonate sferice și cilindrice menționate mai devreme (vezi pagina ) Pentru a determina această metodă, începem prin a sublinia derivarea df/dsf Primul nostru pas este să luăm în considerare o funcție scalară de trei coordonate carteziene/(x, y, z) și să folosim seria Taylor pentru a determina schimbarea în / cauzată de o deplasare de la punctul (x, y, z) la o secundă punctul (x + Дх, у + Ay, z + Az) Găsim pentru această schimbare A/=f A-r+|Av + f Аг + "-' Calculul prezentat aici se referă la o funcție a trei variabile și este o simplă generalizare a calculului de la paginile - , care tratează o funcție a două variabile Gradientul În continuare scriem Д/ în termenii lui As, deplasarea vectorială de la (x, y, z) la (x + Дх, у + Ду, z + Дг) În mod clar (vezi Figura ГѴ- ), Де = i Дх + j Ду + к Дг, astfel încât df df , df\ д + j -f к - j • As + • • • Эх J dy dz J Figura IV- În final, scriem As = u A?, împărțim la As și luăm limita: hm - = - ca ds Mărimea care este punctată în u în această ultimă expresie este apoi recunoscută ca gradient de/în coordonate carteziene Pentru a obține gradientul unei funcții scalare în coordonate cilindrice procedăm aproape în același mod: Considerăm o funcție scalară de trei coordonate cilindrice,/(r, Ѳ, г) Folosind seria Taylor, găsim schimbarea în/ datorită unei deplasări de la punctul (г, Ѳ, z) la un al doilea punct (r + Дг, Ѳ + ДѲ, z + Дг): Д/= у-Дг +ДѲ + - Дг + • • ■ sau ЭѲ dz În continuare, scriem Д/іп termeni ai As Acesta este miezul calculului Din Figura IV- avem As = ёг Дг + ёѳг ДѲ + ё, Дг Gradientul în Cilindrică și sferice Coordonatele Figura IV- Două trăsături ale acestei expresii necesită o discuție În primul rând, deplasarea în direcția creșterii Ѳ (de mărime г ДѲ) este mai degrabă un arc de cerc decât un segment de linie dreaptă Cu toate acestea, din moment ce în cele din urmă vom trece la limită ca As -" , putem considera ДѲ (precum Ar și Дг) ca arbitrar de mici, caz în care arcul este arbitrar aproape de coarda sa de subtensiune Astfel, așa cum este indicat în Figura ГѴ- , Дг, г ДѲ și Дг Figura IV- aproximați la orice grad de precizie dorit trei deplasări reciproc perpendiculare, analogii celor trei deplasări carteziene Ax, Ду și Дг (vezi Figura ГѴ- ) A doua caracteristică a expresiei noastre pentru As care necesită comentariu are de-a face și cu deplasarea în direcția creșterii Ѳ Acesta este: Deoarece arcul face parte dintr-un cerc cu raza r + Дг, ar trebui, strict vorbind, să scriem deplasarea ca (r + Дг) ДѲ, nu г ДѲ Dar termenul suplimentar Дг ДѲ este "de ordinul doi"; adică este produsul a două cantități mici și deci neglijabil în comparație cu г ДѲ Dacă scriem acum expresia noastră pentru termenii Д/іп ai lui As, obținem д/=(ё, |+ё # + ѳ r ЭѲ л у (r)гЭг • Де + ■ • ■ Gradient Notați factorul /r în al doilea termen pentru a compensa factorul r în eer ДѲ în As În cele din urmă, punând As - u As, găsim A/ df L df idf df\ л âs->o As ds \ dr " r ЭѲ dz) Mărimea din expresia precedentă punctată în u este gradientul de/în coordonate cilindrice O procedură analogă poate fi utilizată pentru a găsi gradientul în coordonate sferice; aceasta a fost lăsată ca exercițiu (vezi Problema IV- ) PROBLEME IV- (a) Calculați F = V/ pentru fiecare dintre următoarele funcții scalare: (i) /= xyz (ii) /= x + y + z (iii) f= xy + yz + xz (iv) f= x - z (v) / = e^'siny (b) Verificați acest lucru F • t ds = pentru una sau mai multe dintre funcțiile F determinate în partea (a) alegând pentru curba C: (i) pătratul din planul xy cu vârfurile la ( , ), ( , ), ( , ) și ( , ) (ii) triunghiul din planul yz cu vârfurile la ( , ), ( , ) și ( , ) (iii) cercul cu raza unitară centrat la origine și situat în planul xz (c) Verificați prin calcul direct că VXF = pentru una sau mai multe dintre funcțiile F determinate în partea (a) ГѴ- Verificați următoarele identități în care/și g sunt funcții scalare de poziție diferențiabile arbitrare, iar F și G sunt funcții vectoriale de poziție diferențiabile arbitrare (a) V(/g) = fVg + g Vf (b) V(F • G) = (G • V)F + (F • V)G + FX (VXG) + GX (VXF) (c) V • (/F) =/V • F + F • V/ (d) V • (FXG) = G • (VXF) - F • (VXG) (e) VX (/F) =/VXF + (V/) X F (f) VX (FXG) = (G • V)F - (F • V)G + F(V • G) -G(V • F) (g) VX (VXF) = V(V"F) - V F Probleme IV- Să se arate că VXV/= unde/(x, y, z) este o funcție scalară diferențiabilă arbitrară Să presupunem că derivatele parțiale mixte de ordinul doi sunt independente de ordinea diferențierii De exemplu, d f/dx dz = d fldz dx IV^l (a) Fiecare dintre următoarele funcții este netedă într-o regiune simplu conectată Determinați care dintre ele poate fi scris ca gradient al unei funcții scalare și, pentru cei care pot, utilizați ecuația (IV- ) pentru a găsi acea funcție scalară (i) F = iy- (ii) F = Ck, C o constantă (iii) F = iyz + jxz + kxy (iv) F = ix + jy + kz (v) F = ie~ sin у + je v sin z + kc' sin y (b) Niciuna dintre următoarele funcții nu este netedă peste tot Cu toate acestea, fiecare poate fi scris ca gradient al unei funcții scalare Utilizați ecuația (IV- ) pentru a găsi acea funcție scalară (i) F = r/r , r = ix + jy (ii) F = r/r , r = ix + jy + kz IV- Funcția F(r, Ѳ, z) definită în Problema III- este netedă și are curl zero într-o regiune neconexă simplă constând din tot spațiul tridimensional cu axa z îndepărtată Arătați că nu există o funcție scalară ф astfel încât F = Ѵф evaluând integrala dreaptă a lui F ♦ t de la punctul P,( , - , ) la punctul P ( , , ) pe două căi diferite : CR, partea dreaptă a cercului cu raza situată în planul xy și centrată la origine (vezi figura), și CL, partea stângă a aceluiași cerc Orientați căile așa cum se arată De ce faptul că cele două căi dau rezultate diferite implică că nu există nicio funcție scalară ф astfel încât F = Ѵф? IV- (a) Un dipol electric de putere p situat la origine și orientat în direcția z pozitivă dă naștere unui câmp electrostatic ip E(r, Ѳ, ф) = - ( e, cos ф + e* sin ф) ire r Gradientul Utilizați ecuația (IV- ) pentru a arăta că potențialul dipolului este dat de Ф(г, Ѳ, ф) p cos ф iten r Informații utile: În coordonate sferice, л dr л d m = е^+е*''* + е'Г Пф^ (b) Calculați fluxul câmpului dipol printr-o sferă de rază R centrată la origine (c) Care este fluxul câmpului dipol pe orice suprafață închisă care nu trece prin origine? IV- Iată o "dovadă" că nu există magnetism Una dintre ecuațiile lui Maxwell ne spune asta V • В = , unde В este orice câmp magnetic Apoi, folosind teorema divergenței, găsim Дв-йл Щѵ-в^ о Deoarece В are divergență zero, știm (vezi problema III- ) că există o funcție vectorială, caii it A, astfel încât В = VX A Combinând aceste ultime două ecuații, obținem I J" n • VXA dS = Apoi aplicăm teorema lui Stokes și rezultatul precedent pentru a găsi A • t ds = JJ ii • VXA = Astfel am arătat că circulația lui A este independentă de cale Rezultă că putem scrie A = Ѵф, unde ф este o funcție scalară Probleme Deoarece curba gradientului unei funcții este zero, ajungem la faptul remarcabil că В = VX Ѵф = ; adică toate câmpurile magnetice sunt zero! Unde am gresit? [Preluat de la G Arfken, Amer J Phys , , ( ) ] IV- Legea lui Fick afirmă că în anumite procese de difuzie densitatea de curent J este proporțională cu negativul gradientului densității p; adică J = - к Vp, unde к este o constantă pozitivă Dacă o substanță cu densitatea p(x, y, z, t) și viteza v(x, y, z, t) difuzează conform legii lui Fick, arătați că fluxul este irrotațional (adică VX v = ) IV- (a) O substanță difuzează conform legii lui Fick (vezi Problema IV- ) Presupunând că materia difuză este conservată, deduceți ecuația de difuzie (b) Bacteriile cu densitatea p difuzează într-un mediu conform legii lui Fick și se reproduc cu o rată Xp pe unitate de volum (X este o constantă pozitivă) Arata asta ^ = *V p + Xp IV- (a) Se spune că un fluid este incompresibil dacă densitatea sa p este o constantă (adică este independentă de x, y, z și t) Utilizați ecuația de continuitate pentru a arăta că viteza v a unui fluid incompresibil satisface ecuația V • v = (b) Dacă VX v = , se spune că fluxul de fluid este irrotațional Arătați că pentru un fluid incompresibil care trece printr-un flux irrotațional, Ѵ ф = , unde ф, o funcție scalară numită potențial de viteză, este astfel definită încât v = Ѵф IV- Căldura Q într-un corp de volum V este dată de unde c este o constantă numită căldură specifică a corpului, iar T(x, y, z, t) și p(x, y, z) sunt, respectiv, temperatura și densitatea corpului (Rețineți că presupunem că densitatea este independentă de timp ) Rata cu care căldura curge prin , suprafața de delimitare a corpului, este dată de Gradientul unde к (constantă presupusă) este conductivitatea termică a corpului, iar integrala este preluată pe suprafața care limitează corpul Utilizați aceste fapte pentru a deriva ecuația fluxului de căldură unde a = cp/k IV- În mecanica cuantică non-relatistă, o parte de masă m care se mișcă într-o potențială V(x, y, z) este descrisă de ecuația Schrodinger h дф ~-Ѵ ф + Ѵф = m * ЭГ unde h este constanta lui Planck împărțită la tt și ф(х, у, z, f), care este complexă, se numește funcție de undă Mărimea p = ф*ф este interpretată ca densitate de probabilitate (a) Folosiți ecuația Schrodinger pentru a obține o ecuație de forma Эр ^ + VJ = și obțineți astfel o expresie pentru J în termeni de ф, ф*, m și h (b) Oferiți o interpretare a lui J și a ecuației derivate în (a) IV- (a) Aflați densitatea de sarcină p(x, y, z) care produce câmpul electric E = g(ix + jy + kz), unde g este o constantă (b) Găsiți o potențială electrostatică Ф astfel încât - ѴФ este câmpul E dat în (a) (c) Verificați că Ѵ Ф = - p/e IV- (a) Pornind de la teorema divergenței, derivați ecuația lui V u + (Vu) ■ (Vt>)] dV, unde и și v sunt funcții scalare de poziție și S este o suprafață închisă care cuprinde volumul V Aceasta este uneori numită prima formă a teoremei lui Green (b) Dacă V u = folosiți prima formă a teoremei lui Green pentru a arăta că fi n ■ (u Vu) dS - fflrlv Nv unde |Ѵи| = (Vu) • (Ѵи) (c) Folosiți prima formă a teoremei lui Green pentru a arăta că s v Ѵи) dS = JJJ (и V u - v Ѵ и) dV Aceasta este a doua formă a teoremei lui Green Probleme IV- O ecuație de formă ѵу- Й V dt unde f este o funcție de două ori diferențiabilă a poziției și timpului, se numește ecuație de undă Descrie o undă care se propagă în spațiu cu viteza v Folosiți ecuațiile lui Maxwell (Problema III- ) pentru a arăta că, în absența sarcinilor și a curenților (adică, p și J ambele zero), toate cele trei componente carteziene ale lui E și В satisface o ecuație de undă cu v = c, unde c = /Veop,o este viteza luminii De exemplu, Ѵ ЕД = d Ex гіё' Astfel, existența undelor electromagnetice care călătoresc în spațiul gol cu viteza luminii este o consecință a ecuațiilor lui Maxwell IV- (a) În text am găsit potențialul și câmpul pentru cazul unui cilindru infinit între plăci paralele cu cilindrul menținut la potențialul zero Cum trebuie modificată soluția dacă cilindrul este ținut la o potențială Vo #= ? (b) Arătați că nu există nicio sarcină netă pe cilindru IV- (a) O sferă cu raza R este situată între două plăci paralele foarte mari care sunt separate de o distanță s Se menține o diferență de potențial între plăci, iar sfera este menținută la potențialul zero Găsiți potențialul și câmpul peste tot în afara sferei și între plăci Să presupunem că R" s (b) Arătați că nu există nicio sarcină netă pe sferă (c) Repetați partea (a) presupunând că sfera este ținută la o potențială Vo*O IV- Fie f(x, y) o funcție scalară diferențiabilă a lui x și y și fie u = i cos Ѳ + j sin Ѳ Transformați într-un sistem de coordonate rotit x', y' astfel încât x' să fie paralel cu u (vezi figura) Arătați că derivata direcțională în direcția lui u este dată de df - vr % ds dx Gradientul IV- Sunteți într-un punct (a, b, c) de la suprafață z = (r - x - y )'/ (z > ) Presupunând că ambele a și b sunt pozitive, în ce direcție trebuie să vă deplasați (a) astfel încât rata de schimbare a lui z să fie zero? (b) astfel încât rata de creștere a lui z să fie cea mai mare? (c) astfel încât rata de scădere a lui z să fie cea mai mare? Desenați o schiță pentru a arăta semnificația geometrică a răspunsurilor dvs ГѴ- Vectorul unitar normal la suprafata z = fix, y) este dat de [vezi ecuația (П- )] De asemenea, am stabilit că VF este un vector normal la suprafață Fix, y, z) = const (pagina ) astfel încât VF/|VF| este un vector unitar normal la suprafață Fix, y, z) = const Să se arate că aceste două expresii pentru vectorul normal unitar sunt identice dacă F(x,y,z) = const și z = fix, y) descriu aceeași suprafață IV- Utilizați rezultatele problemei - și expresia gradientului în coordonate cilindrice (vezi pagina ) pentru a obține forma laplacianului în coordonate cilindrice, dată la pagina IV- Utilizând procedura descrisă în text (paginile - ) obțineți expresia pentru gradientul lui ф în coordonate sferice: л дф л дф л Эф Ѵф-е, -+еФ -+еѳ - ІѴ- Folosiți rezultatele problemei - și expresia gradientului în coordonate sferice derivată în problema IV- pentru a obține forma laplacianului în coordonate sferice dată la pagina IV- Să presupunem că găsiți o soluție a ecuației lui Laplace care satisface anumite condiții la limită Este această soluție unică sau există și altele? Această problemă va răspunde la această întrebare în anumite cazuri simple Se consideră regiunea spațiului complet închisă de o suprafață So și care conține în interiorul ei obiectele , , , (dintre care două sunt ilustrate în diagramă) Să presupunem că So se menține la o potențială constantă Фо, obiectul nr la Ф], obiectul nr la Ф și așa mai departe Apoi, în regiunea fără încărcare F închisă de So și între obiecte, potențialul trebuie să satisfacă ecuația lui Laplace Ѵ Ф = și condițiile de limită Ф = Фо pe So Ф, pe Sj Ф pe S Probleme Următorii pași vă vor ghida printr-o dovadă că Ф este unică (a) Să presupunem că există două potențiale и și v, ambele satisfăcând ecuația lui Laplace și condițiile la limită enumerate mai devreme Formați diferența lor w = и - v Arătați că V w = în R (b) Care sunt condițiile la limită îndeplinite de w? (c) Aplicați teorema divergenței la n • (w Vw) dS, unde integrarea se realizează pe suprafaţa So + , + S + • • şi se arată prin aceasta că dV = , unde V este volumul regiunii R (d) Din rezultatul (c) argumentați că Vw = și că aceasta, la rândul său, înseamnă că w este o constantă (e) Dacă w este o constantă, care este valoarea acesteia? (Folosiți condițiile la limită de pe w pentru a răspunde la aceasta ) Ce spune aceasta despre и și u? (f) Dovada unicității prezentată la (a) la (e) implică specificarea valorii potențialului pe diferite suprafețe Am fi putut specifica un alt tip de condiție de limită și am fi dovedit totuși unicitatea? Dacă da, în ce fel sau moduri ar diferi dovada și rezultatul de cele date mai sus? IV- În text am definit gradientul în termenii anumitor derivate parțiale Este posibil să oferim o definiție alternativă similară ca formă cu definițiile noastre ale divergenței și buclei Prin urmare, V/= lim - - [ f îifdS J дѵ-"о AV J Js Aici/este o funcție scalară a poziției, S o suprafață închisă și AU volumul pe care îl înconjoară Ca de obicei, n este un vector unitar normal la și care arată din volumul inclus (a) Urmând o procedură similară cu cea folosită în text în tratarea divergenței, integrați peste un "cuboid" și arătați că definiția anterioară dă expresia df df df V/=i^ +j^ + k/ J dx J dy dz Gradientul (b) Utilizați definiția alternativă a gradientului dată mai sus pentru a arăta că derivata direcțională a/în direcția specificată de vectorul unitar u este dată de df - w - = u • Nf ds [Sugestie: Evaluează и • j js *fdS = jj " " peste un cilindru mic (lungime As, aria secțiunii transversale ДА; vezi figura) a cărui axă este în direcția vectorului unitar constant u Apoi împărțiți la volumul cilindrului (Ca ДА) și luați limita pe măsură ce volumul se apropie de zero ] (c) Argumentând așa cum am făcut în text pentru a stabili teorema divergenței, utilizați definiția alternativă a gradientului pentru a arăta că unde este o suprafață închisă care cuprinde volumul V (d) Obţineţi relaţia enunţată la (c) direct din teorema divergenţei [Sugestie: În JfsF • n dS = fffvV • F dVput F = e/ unde e este un vector unitar constant ] Probleme (e) Verificați relația menționată în (c) pentru funcția scalară /=x +y + z integrându-se peste cilindrul unitar prezentat în figură ГѴ- (a) Se consideră o suprafață z = f(x, y) Fie u un vector de lungime arbitrară tangent la suprafață într-un punct P(x, y, z) în direcția vectorului unitar p = ip, + jpy așa cum este indicat în figură Folosiți derivata directă pentru a arăta asta u = p + k(p • V/), unde V/' este evaluat la (x, y) [Notă: Deoarece lungimea lui u este arbitrară, rezultatul dvs poate diferi de cel precedent printr-o constantă multiplicativă pozitivă ] (b) Fie v un al doilea vector de lungime arbitrară tangent la suprafața de la P, dar în direcția vectorului unitar q = \qx + j?v (p + q) Atunci de la (a) avem v = q + k(q"V/) Arata asta u X v = [k • (p X q)](k - V/) și utilizați aceasta pentru a rederive ecuația (П- ) pentru vectorul unitar n normal la suprafața z = /(x, y) la (x, y, z) Aceasta arată că rezultatul derivat în text pentru n este unic (în afară de semn) chiar dacă a fost obținut cu opțiunile speciale p = i și q = j ZV- (a) Folosind ecuațiile lui Maxwell (vezi problema III- ), arătați că putem scrie В = VXA, E = -ѴФ - m unde A (numit potențialul vectorial) este o funcție vectorială a poziției și timpului, iar Ф (potențialul scalar) este o funcție scalară Gradientul funcţie de poziţie şi timp, cu condiţia ca A şi Ф să satisfacă ecuaţiile Ѵ?Ф + A) = p/e°' V A - poeo = -p,oJ + VV • A + p,oeo ~ (b) Arătați că dacă deținem două potențiale noi A' = A + Vx, ф' = ф-^, dt unde x este o funcție scalară arbitrară a poziției și timpului, atunci В = VX A', E = -ѴФ' dt Adică câmpurile E și В nu sunt modificate de modificarea potențialelor A și Ф Modificarea de la (А, Ф) la (А', Ф') se numeşte transformare gauge (c) Să se arate că dacă avem nevoie de x pentru a satisface ecuația V"A + eop,o - , d x VX - eoM-o = - apoi Г л' ЭФ' nv • A + e |x = °- (d) Dacă x satisface ecuația dată în (c), arătați că A' și Ф' satisfac ecuațiile e și , A' t VA - eo|xo-=-p-oJ Ideea pentru toate acestea este că putem face o transformare gauge [ca în (b)], impunem condiția dată în (c) și, prin urmare, obținem o potențială scalară și vectorială care satisface ecuațiile din partea (d), care sunt ecuații de undă cu termeni sursă proporționali cu p și J TV- Se poate scrie ecuația de mișcare a unui fluid ideal + (v • V)v dV ot unde V este volumul fluidului și S este suprafața acestuia Неге feM(x, y, z) este forța externă pe unitatea de masă care acționează asupra fluidului, p(x, y, z) este Probleme presiunea fluidului și p(x, y, z) este densitatea acestuia, toate într-un punct (x, y, z) al fluidului, iar v(x, y, z, t) este viteza fluidului în punctul (x, y, z) și la momentul t (a) Folosiți forma teoremei divergenței dată în problema IV- (c) pentru a rescrie ecuația de mișcare a unui fluid ideal sub forma fext-|Vp = |+(vV)v (b) Să se arate că în cazul static (v = ), ecuația mișcării devine fext = ( /p) Vp (c) Luați în considerare o coloană de fluid incompresibil orientată vertical paralelă cu axa z așa cum se arată în figură Presupunând că singura forță externă care acționează asupra fluidului este atracția gravitațională uniformă descendentă a pământului, aplicați ecuația pentru cazul static din (b) pentru a arăta că P = Po " Pgz unde g este accelerația datorată gravitației și p este o constantă A Soluții la probleme Trebuie să înveți făcând lucrul; căci, deși crezi că știi, nu ai nicio certitudine până nu încerci Sofocle Capitolul I (a) (ix + jy)/Vx + у (b) (i + j)(x + y) /V (c) -iy + jx (d) (ix + jy + kzj/Vx + у + г • (a) (a cos ыі + b sin /> cos cot (viteza) -ico a cos (b) E = roi/eo, x Î și - roi/eo, x ȘI (Cu) (A) (b) (c) (A) (c) irR e~R (d) [E(b) - E( )]d și -cri/ e , x хі/е , - b SxSb; )i (+ forx > , - forx b și (p r/ e )ez , r b, (Po/ eo)ro, r Ь Câmpul este zero pentru r > b dacă p, = - p / Încărcarea totală este atunci zero Е = (с) Е = Rezolvarea problemelor (a) (x + у + z) (e) -y/tf + у ) (b) (f) (c) -(e"x + e~y + е г) (g) (d) z (h) Integrala de suprafață este egală cu s (x + y + zo) pentru funcția problemei II- a Integrala de suprafață este egală cu pentru funcția problemei П- Ь Integrala de suprafață este egală cu s (e~sri - es/ )(e-^ + е~Уо + е г") pentru funcția (II- c) (b) V • G = f(f ) = constant/r (a) ЗВ (b) tt/? ?/ (c) irB (b) ітЛ В Capitolul III (a) (-iy + jz + kx) (e) -ix - jy + kz (b) jx (f) (iy - jx) (c) ie~z + je~x + key~y (g) iz - kx (d) (h) Integrala dreaptă este egală cu x v pentru funcția Problemei III- a Integrala dreaptă este egală cu pentru funcția Problemei III- b Integrala dreaptă este egală cu s(esl - e~s' )e~y° pentru funcția problemei Ш-Зс Integrala dreaptă este egală cu pentru funcția problemei ПІ- d (a) a / + a / / l a | l (b) г -> + a = as а V \ / Ѵз (c) n • VXF = (l + y) = -at у = ѴЗ V (d)and( h) (a) Integrală de linie și integrală de suprafață egală cu (b) Integrală de linie și integrală de suprafață egală -Зтг/ (c) Integrală de linie și integrală de suprafață egală cu - irB / x (c) V • G (d) H = -j(x - x )z + k(x + x )y (e) V • G + [Notă: rezultatele dvs pot diferi de acestea prin constante aditive ] Solutii la Probleme (a) $(H • t ds = ffsG • n dS unde S este o suprafață de acoperire a curbei închise C (d) Suprafața și integrala de volum sunt egale fiecare cu i - j - k Capitolul IV (a) (i) F = iyz + jxz + kxy (ii) F = (Lc + jy + kz) (iii) F = i(y + z) + j(x + z) + k(x + y) (iv) F = ix - kz (v) F = -ie~x sin у + cos y (a) (i) Nu este independent de cale (ii) ф = cz + const (iii) ф = xyz + const (iv) ф - |(x + y + z ) + const (v) Nu este independent de cale (b) (i) ф = In r + const (ii) ф = |r / + const (a) p = ge (b) Ф = -^(x + y + z ) (a) Adăugați Vo la rezultatul obținut în text (а) Ф(г ф) = - Eor(l - R'lri) cos ф unde sfera este centrată la origine şi cele două plăci sunt paralele cu planul xy şi situate la z - ±s! (c) Adăugați Vo la rezultatul dat în partea (a) (a) Deplasați-vă în direcția ±(i/> - ja)/Va + b (b) Deplasați-vă în direcția gradientului: -(ia + jZ>)/Vr - a - b (c) Deplasați-vă în direcția opusă gradientului: (ia + j/?)/Vr - a - b (e) Integrala de suprafață și integrala de volum sunt egale fiecare cu тгк Index Legea circuitului lui Ampere, - formă diferențială a, Arfken, G" Amold, Matei, Unghi azimutal, ix Beethoven, Ludwig van, v Suprafata de acoperire, - și teorema lui Stokes, - Forța centrală, și funcția iritațională Cervantes, Miguel de, citat, ix Chalmers, PR, Densitatea sarcinii, Circulație, și curl, Legea circulației, formă diferențială a, - Conservare, - Ecuația de continuitate, Forța Coulomb, , Legea lui Coulomb, - și independența căii, - "Cuboid", n cilindric, ^ Curl, - , definiție alternativă a, - în coordonate carteziene, și circulație, în coordonate cilindrice și sferice, - sensul lui, - ca operator, și independența căii, - Densitatea curentului, Cuboid cilindric, - Del notation, - , - curl operator in, divergence operation in, gradient operation in, and Laplacian, - Denham, Sir John, citat, Density, charge, curent, Ecuații diferențiale, n Ecuație de difuzie, Curbă direcționată, Direcțional derivat, - Divergence, ^ în coordonate carteziene, - , în coordonate cilindrice și sferice, - , și operator, Teorema divergenței, - , , aplicații ale, - derivație - ilustrare a, ^ declarație a, Teorema Stones, reacția la, - în două dimensiuni, valabilitate - de, - Sarcina electrică, - , Câmp electric, ca funcție vectorială, Teoria electromagnetică și ecuațiile lui Maxwell, Index Forța electromotoare, Câmp electrostatic, determinarea - a, - și legea lui Gauss, - și integrală de linie, - și funcție vectorială, Potențialul electrostatic, Electrostatică, - și calcul vectorial, - Ecuația mișcării, a fluidului ideal, - Faraday, Michael, Legea lui Fick, Câmp, electrostatic, - Linii de câmp, Teoria câmpului, electromagnetism ca, Fluid, ideal, - incompresibil, Ideal fluid, ecuația - Flux, - , - "Formal", Funcția n din mai mult de o variabilă, dintr-o variabilă, - Transformarea ecartamentului, Legea lui Gauss, - , și suprafață închisă, formă diferențială a, , , , - și teorema divergenței, utilizare în găsirea câmpului, - și flux, - teorema lui Gauss Vezi teorema divergenței Semnificația geometrică a gradientului, - Gradient, - definiție alternativă a, în coordonate carteziene, în coordonate cilindrice și sferice, - și derivată direcțională, - în găsirea câmpului electrostatic, - semnificația geometrică a, - și ecuația lui Laplace, - ca operator, Teorema lui Green, - Căldură, Fluxul de căldură, analiza, ecuația , - Fluid ideal, ecuația de mișcare a, - Flux iritațional, In, Suprafata izoterma, Ecuația lui Laplace, - în coordonate carteziene, în coordonate cilindrice, Soluții ale, - în coordonate sferice, unicitate ale Soluțiilor, , - Laplacian, - ca operator, Integrale drepte, - definiția, - evaluarea, - , - care implică funcția vectorială, - și independența de cale, - , - Camp magnetic și legea circuitală a lui Ampere, - dovada absenței, - Matematică și contexte non-matematice, Ecuațiile lui Maxwell, și undele electromagnetice, și legea lui Gauss, și teorema lui Stokes Unități MKS, Mișcarea fluidului ideal, - Vector normal, - formă alternativă a, , unicitate a, , unitate normală, Suprafata deschisa, Operator, Condensator cu plăci paralele, - Independența traseului, - și forțele centrale, - și curl, - și integrală de linie, - , - Index Permitivitatea spațiului liber, constanta lui Planck, ecuația lui Poisson, Unghiul polar, ix Principiul suprapunerii, , Vector radial, Unități MKS raționalizate, regulă pentru mâna dreaptă, - "Rotație", in Funcția scalară, Potențialul scalar, - Ecuația Schrodinger, regiuni pur și simplu conectate și Teorema lui Stokes, - , - Funcția Smooth, n, Funcția vectorială solenoidală, Sofocle, citat, , Condensator sferic, - Teorema lui Stokes, - aplicații ale, - derivarea, - teorema de divergență, relația cu, , - ilustrare a, - și regiunile pur și simplu conectate, - , - în două dimensiuni, valabilitatea, Indice, n Suprafețe de subvolum, Suprapoziție, principiul, , Integrală de suprafață, definiția, - evaluarea, - , Vector tangent, - seria Taylor, - Torus, Vector normal unitar Vezi vectorul normal Calcul vectorial și electrostatică, Funcții vectoriale, - și câmp electrostatic, de gradienți ca, integrale de linie care implică, - de poziție, solenoidal, Potențial vectorial, Câmp de viteză - Potențial de viteză, Ecuația de undă, Funcția de undă, Lucrare, , - CARTESIAN z CILINDRIC DIVERGENȚĂ div FV*F GRADIENT grad/ V/ dFx dFy dF, ore dx dy dz până la dFa dF, -fTr(rF'> + f"M+^Z = eu - Tr CURL curl FVXF LAPLACIAN V / (Vf) = t'L VJ ,y dy УУ) = VJ)z dz m = î (V (V (V d Iy dF dy dz dFx dF F^y~ dz ~dx F>-^' F) (V (V dx dy (V F), F)e F> dx dy dz Â Г dr -idFț-dF, г дѲ dz dh dr dF, dz in absenta Г dr dFr (rF^) dr d f d / Г дѲ dz SFERIC dF г sin ЭѲ df TM = /г df df (^f )er sjn ф (VX Е)ф dFit (VXF), = -sjn~j~ ^ф (sm Ф Faith) г Sin ЭѲ - X Э z р х г sin ЭѲ г дг ѳ până la dFr = ~г~дг " aug (VX В )ѳ д r dg dr -^- ^f + G sin f df \ °F / df\ , d f Eph ) r sin f dѲ BIBLIOTECA PUBLICA CHICAGO ROSO CCHD HWLC ETAJ 